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Abstract. We present an overview of recent results and developments of the no-core 
shell model (NCSM), an ab initio approach to the nuclear many-body problem for 
light nuclei. In this approach, we start from realistic two-nucleon or two- plus three- 
nucleon interactions. Many-body calculations are performed using a finite harmonic- 
oscillator (HO) basis. To facilitate convergence for realistic inter-nucleon interactions 
that generate strong short-range correlations, we derive effective interactions by unitary 
transformations that are tailored to the HO basis truncation. For soft realistic 
interactions this might not be necessary. If that is the case, the NCSM calculations 
are variational. In either case, the ab initio NCSM preserves translational invariance 
of the nuclear many-body problem. In this review, we, in particular, highlight results 
obtained with the chiral two- plus three-nucleon interactions. We discuss efforts to 
extend the applicability of the NCSM to heavier nuclei and larger model spaces using 
importance-truncation schemes and/or use of effective interactions with a core. We 
outline an extension of the ab initio NCSM to the description of nuclear reactions by 
the resonating group method technique. A future direction of the approach, the ab 
initio NCSM with continuum, which will provide a complete description of nuclei as 
open systems with coupling of bound and continuum states, is given in the concluding 
part of the review. 
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1. Introduction 

A major outstanding problem in nuclear physics is to calculate properties of finite nuclei 
starting from the basic interactions among nucleons. This problem has two parts. First, 
the basic interactions among nucleons are complicated. They are not uniquely defined 
and there is clear evidence that more than just two-nucleon forces are important. Second, 
the nuclear many-body problem is very difficult to solve. This is a direct consequence 
of the complex nature of the inter-nucleon interactions. Both short-range and medium- 
range correlations among nucleons are important and for some observables long-range 
correlations also play a significant role. 

In this review, we focus mainly on the second part of the problem, namely on 
the solution of the many-nucleon problem. The two-nucleon interactions we take as 
an input provided to us by other theorists. We do, however, utilize three- and many- 
nucleon calculations to determine parameters of three-nucleon interactions. We also 
note that, in general, few-nucleon and many-nucleon calculations provide feedback to 
those constructing potentials. 

Various methods have been used to solve the few-nucleon problem in the past. 
The Faddeev method [1] has been successfully applied to solve the three-nucleon 
bound-state as well as the scattering problem for different nucleon-nucleon (NN) 
potentials [21 [3l H]. For the solution of the four-nucleon problem one can employ 
Yakubovsky's generalization of the Faddeev formalism [5], as done, e.g., in Refs. [6] or 
[7j. Alternatively, other methods have also been succesfuUy used, such as, the correlated 
hyperspherical harmonics expansion method [8l [9] or the Green's function Monte Carlo 
method (GFMC) [10]. Recently, a benchmark calculation by seven different methods 
was performed for a four-nucleon bound state problem [TT] giving the same result within 
error. However, there are few approaches that can be successfully applied to solve the 
bound-state problem in systems of more than four nucleons, when realistic inter-nucleon 
interactions are used. These include the Green's function Monte Carlo method, which is 
capable of solving the nuclear many-body problem with realistic interactions for systems 
of up to A = 12 and the coupled cluster method [121 [131 [HI [13 [IS] , which is applicable 
typically to closed-shell and nearby nuclei. 

The solution of the nuclear many-body problem is still more complex when 
scattering or nuclear reactions are considered. For A = 3 and 4 nucleon systems, the 
Faddeev and Faddeev- Yakubovsky as well as the hyperspherical harmonics (HH) [E] 
or the Alt, Grassberger and Sandhas (AGS) [H] methods are applicable and successful. 
However, ab initio calculations for scattering processes involving more than four nucleons 
overall are challenging and still a rare exception |19j . 

Nuclei are open systems with bound states, weakly bound halo states, unbound 
resonances as well as scattering states. A realistic ab initio description of light nuclei 
with predictive power must have a capability to describe all the above classes of states 
within a unified framework. Coupling to the continuum cannot be neglected. 

In this review, we describe the ab initio no-core shell model (NCSM) [20], another 
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method for solving the nuclear many-body problem developed recently and applicable 
to light nuclei up to A = 16 and beyond. The first no-core shell model calculations [2T| 
were performed with G- matrix-based two-body interactions [22] . Later, the Okubo-Lee- 
Suzuki procedure [23l [2l] was implemented to derive two-body effective interactions 
for the NCSM [25]. This resulted in the elimination of the purely phenomenological 
parameter used to define the G-matrix starting energy. A truly ab initio formulation 
of the approach was presented in Ref. [26], where convergence to the exact bound-state 
solutions was demonstrated for the A = 3 system. Here, in the first part of the review, 
we discuss the NCSM in its standard formulation, applicable to the nuclear bound- 
state problem. In Sect. O we briefly present the NCSM formalism. In Sect. El we 
show recent results obtained with the chiral NN plus three-nucleon (NNN) interactions. 
Calculations of radii, moments and transitions of He, Li and Be isotopes using different 
realistic NN potentials are discussed in Sect. HI In Sect. [5], we describe efforts to extend 
the applicability of the NCSM to heavier nuclei. In the second part of the review in 
Sections [6HHI we discuss in detail new developments and outline future efforts to extend 
the NCSM by including continuum states to describe unbound states, scattering and 
nuclear reactions in a unified framework. Conclusions are given in Sect. M 

2. Ab initio no-core shell model 

In the ab initio no-core shell model, we consider a system of A point-like non-relativistic 
nucleons that interact by realistic two- or two- plus three-nucleon interactions. By 
the term "realistic two-nucleon interactions", we mean NN potentials that fit nucleon- 
nucleon phase shifts with high precision up to a certain energy, typically up to 350 
MeV. A realistic NNN interaction includes terms related to two-pion exchanges with 
an intermediate delta excitation. In the NCSM, all the nucleons are considered active, 
there is no inert core like in standard shell model calculations. Hence, the "no-core" in 
the name of the approach. 

There are two other major features in addition to the employment of realistic 
NN or NN-I-NNN interactions. The first one is the use of the harmonic oscillator 
(HO) basis, truncated by a chosen maximal total HO energy of the A-nucleon system. 
The reason behind the choice of the HO basis is the fact that this is the only basis 
that allows for the use of single-nucleon coordinates and, consequently, the second- 
quantization representation, without violating the translational invariance of the system. 
The powerful techniques based on the second quantization and developed for standard 
shell model calculations can then be utilized. Therefore, the "shell model" in the name 
of the approach. As a downside, one has to face the consequences of the incorrect 
asymptotic behavior of the HO basis. 

The second feature comes as a result of the basis truncation. Standard, accurate 
NN potentials, such as the Argonne V18 (AV18) [27], CD-Bonn 2000 [28], INOY (inside 
non-local outside Yukawa) [29] and, to some extent, also the chiral N^LO [30], generate 
strong short-range correlations that cannot be accomodated even in a reasonably large 
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HO basis. In order to account for these short-range correlations and to speed up 
convergence with the basis enlargement, we construct an effective interaction from the 
original, realistic NN or NN+NNN potentials by means of a unitary transformation. 
The effective interaction depends on the basis truncation and by construction becomes 
the original, realistic NN or NN+NNN interaction as the size of the basis approaches 
infinity. 

Recently, a new class of soft potentials has been developed, mostly by means of 
unitary transformations of the standard, accurate NN potentials mentioned above. 
These include the Viowk [M], the Similarity Renormalization Group (SRG) [32] and 
the UCOM ^3] NN potentials. A different class of soft phenomenological NN potential 
used in some NCSM calculations are the JISP potentials based on inverse scattering [3^ . 
These soft potentials are to some extent already renormalized for the purpose of 
simplifying many-body calculations. Therefore, we can perform convergent NCSM 
calculations with these potentials unmodified, or "bare." In fact, the chiral N^LO NN 
potential [30] can also be used bare with some success. NCSM calculations with bare 
potentials are variational with the HO frequency and the basis truncation parameter as 
variational parameters. 



2.1. Hamiltonian 



The starting Hamiltonian of the ah initio NCSM is 

^ _ ^ y A A 

^ ]4 5Z ^ VNN,ij + ^ ^NNN.iifc , (1) 

i<j i<j i<j<k 

where m is the nucleon mass, VNN,ii is the NN interaction, and \4fNN,ijfc is the three- 
nucleon interaction. In the NCSM, we employ a large but finite HO basis. When soft 
NN potentials are used, it is often feasible to employ a sufficiently large basis to reach 
convergence with the Hamiltonian ([1]). 

On the other hand, if realistic nuclear interactions that generate strong short-range 
correlations are used in Eq. ([T]), we must derive an effective interaction appropriate 
for the basis truncation. To facilitate the derivation of the effective interaction, we 
modify the Hamiltonian ([T]) by adding to it the center-of-mass (CM) HO Hamiltonian 
HcM = TcM + UcM, where Ucu = lAmVl^R'^, R = \ n. The effect of the HO CM 
Hamiltonian will later be subtracted out in the final many-body calculation. Due to the 
translational invariance of the Hamiltonian ([1]), the HO CM Hamiltonian has in fact no 
effect on the intrinsic properties of the system. The modified Hamiltonian can be cast 
into the form 

A 



Ha = Ha + Hi 
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i<j<k 



(2) 



2.2. Basis 

In the ab initio NCSM, we use a HO basis that allows preservation of translational 
symmetry of the nuclear self-bound system, even if single-nucleon coordinates are 
utilized. This is possible as long as the basis is truncated by a maximal total HO 
energy of the A-nucleon system. A further advantage is that the HO wave functions 
have important transformation properties [35] that facilitate and simplify calculations. 
A single-nucleon HO wave function can be written as 

^nim{r] b) = Rni{r; b)Yira{f) , (3) 

with Rni{r, b), the radial HO wave function, and b, the HO length parameter related to 



the HO frequency as 6 = y with m the nucleon mass. 

Because the NN and NNN interactions depend on relative coordinates and/or 
momenta, the natural coordinates in the nuclear problem are the relative, or Jacobi, 
coordinates. For the present purposes we consider just a single set of Jacobi coordinates 
(a more general discussion can be found in Ref. [36]): 
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(4) 
(5) 
(6) 

(7) 



Here, is proportional to the center of mass of the A-nucleon system. On the other 
hand, ^ is proportional to the relative position of the p + 1-st nucleon and the center 
of mass of the p nucleons. 



2.2.1. Antisymmetrization of Jacobi- coordinate HO basis As nucleons are fermions, 
we need to construct an antisymmetrized basis. The way to do this, when the Jacobi- 
coordinate HO basis is used, is extensively discussed in Refs. [261 EHl [37] . Here we briefly 
illustrate how to do this for the simplest case of three nucleons. 

One starts by introducing a HO basis that depends on the Jacobi coordinates 
and ^2, defined in Eqs. (5) and (6), e.g., 

\{nls]t-MCJ)JT) . (8) 

Here n, / and A/", C are the HO quantum numbers corresponding to the harmonic 
oscillators associated with the coordinates (and the corresponding momenta) ^ and 
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^2, respectively. The quantum numbers s,t,j describe the spin, isospin and angular 
momentum of the relative-coordinate two-nucleon channel of nucleons 1 and 2, while JT" 
is the angular momentum of the third nucleon relative to the center of mass of nucleons 1 
and 2. The J and T are the total angular momentum and the total isospin, respectively. 
Note that the basis (IHl) is antisymmetrized with respect to the exchanges of nucleons 
1 and 2, as the two-nucleon channel quantum numbers are restricted by the condition 
(— 1)'^''^* = — 1. It is not, however, antisymmetrized with respect to the exchanges of 
nucleons 1^3 and 2^3. In order to construct a completely antisymmetrized basis, 
one needs to obtain eigenvectors of the antisymmetrizer 



where T^~^^ and T^~^ are the cyclic and the anti-cyclic permutation operators, 
respectively. The antisymmetrizer A" is a projector satisfying XX = X. When 
diagonalized in the basis ([H]), its eigenvectors span two eigenspaces. One, corresponding 
to the eigenvalue 1, is formed by physical, completely antisymmetrized states and the 
other, corresponding to the eigenvalue 0, is formed by spurious states. There are about 
twice as many spurious states as the physical ones [38] . 

Due to the antisymmetry with respect to the exchanges 1^2, the matrix elements 
in the basis (IHl) of the antisymmetrizer X can be evaluated simply as {X) = |(1 — 2P2,3), 
where P2,3 is the transposition operator corresponding to the exchange of nucleons 2 and 
3. Its matrix element can be evaluated in a straightforward way (see e.g., Ref. [26] ) 



where Ni = 2ni + k + 2Mi + Ci,i = 1,2; j = V2j + 1; and {nihAfiCiL\Af2C2n2l2L)3 
is the general HO bracket for two particles with mass ratio 3, as defined, e.g., in Ref. 
[39]. The expression ffTOj) can be derived by examining the action of P2,3 on the basis 
states dH]). That operator changes the state \nl{^i),J\fC{^2),L) to \nl{^\),J\fC{^'2), L) , 
where = 1,2 are defined as ^i,i = 1,2 but with the single-nucleon indexes 2 
and 3 exchanged. The primed Jacobi coordinates can be expressed as an orthogonal 
transformation of the unprimed ones, see e.g., Ref. |26j. Consequently, the HO wave 
functions depending on the primed Jacobi coordinates can be expressed as an orthogonal 
transformation of the original HO wave functions. Elements of the transformation are 
the generalized HO brackets for two particles with the mass ratio d, with d determined 
from the orthogonal transformation of the coordinates, see e.g. Ref. [35] . 



A' = -(l + T(-) + T(+)) , 



(9) 



{{nihSiJiti-J\flCiJi)JT\P2,3\{n2l2S2]2t2\M2C2j2)JT) 




(10) 
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The resulting antisymmetrized states can be classified and expanded in terms of 
the original basis ([8]) as follows 

\mjT) = ^{nls]t-MCJ\\NiJT)\{nls]t-MCJ)JT) , (11) 

where N = 2n + I + 2M + C and where we have introduced an additional quantum 
number i that distinguishes states with the same set of quantum numbers J,T, e.g., 
z = 1, 2, . . . r with r the total number of antisjTiimetrized states for a given A^, J, T. The 
symbol {nlsjt]AfCJ'\\NiJT) is a coefficient of fractional parentage. 

2. 2. 2. Slater determinant basis A generalization to systems of more than three nucleons 
can be done as shown, e.g. in Ref. [36]. It is obvious, however, that as we increase 
the number of nucleons, the antisymmetrization becomes more and more involved. 
Consequently, in standard shell model calculations one utilizes antisymmetrized wave 
functions constructed in a straightforward way as Slater determinants of single-nucleon 
wave functions depending on single-nucleon coordinates ipi{fi). It follows from the 
transformations of HO wave functions that the use of a Slater determinant basis 
constructed from single nucleon HO wave functions, such as, 

y^nijmmtir, (T, T] h) = Rni{r; b){Yi{f)x{a)YZh{T)mt , (12) 

results in eigenstates of a translationally invariant Hamiltonian that factorize as products 
of a wave function depending on relative coordinates and a wave function depending 
on the CM coordinates. This is true as long as the basis truncation is done by a 
chosen maximum of the sum of all HO excitations, i.e., '^f^ii^Ui + /j) < Ntotmax- In 
Eq. ( JT2l) . a and r are spin and isospin coordinates of the nucleon, respectively. The 
physical eigenstates of a translationally invariant Hamiltonian can then be selected as 
eigenstates with the CM in the state: 

(ri . . . taOi . . . a An . . . ta\AXJMTMt)sd 
= (fi . . . eVi^Ti . . . (TATi . . . TA\AXJMTMT)ipoooi^o; b) . (13) 

For a general single-nucleon wave function this factorization is not possible. The use 
of any other single-nucleon wave function than the HO wave function will result in the 
mixing of CM and internal motion. 

In the ab initio NCSM calculations, we use both the Jacobi-coordinate HO basis 
and the single-nucleon Slater determinant HO basis. One can choose whichever is more 
convenient for the problem to be solved. One can also mix the two types of bases. In 
general, for systems of A < 4, the Jacobi coordinate basis is more efficient, as one can 
perform the antisymmetrization easily. The CM degrees of freedom can be explicitly 
removed and a coupled J^T basis can be utilized with matrix dimensions of the order 
of thousands. For systems with A > A, it is, in general, more efficient to use the Slater 
determinant HO basis. In fact, we use the so-called m-scheme basis with conserved 
quantum numbers M = rrii, parity vr and Mt = X^iLi ''^ti- The antisymmetrization 
is trivial, but the dimensions can be huge, as the CM degrees of freedom are present, and 
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no JT coupling is considered. The advantage is the possibihty to utihze the powerful 
second-quantization technique, shell model codes, transition density codes and so on. 

As mentioned above, the model space truncation is always done using the condition 
X;^i(2ni + k) < Ntotmax- Often, instead of Ntotmax, we introduce the parameter A^^ax 
that measures the maximal allowed HO excitation energy above the unperturbed ground 
state. For A = 3,4 systems Aniax = Ntotmax- For the p-shell nuclei they differ, e.g. for 

^Li, iVmax = Ntotmax ~ 2, for ^^C, A^ax = ^totmax ~ 8) C^C. 

2.3. Effective interaction 

In the ah initio NCSM calculations we use a truncated HO basis, as discussed in previous 
sections. The inter- nucleon interactions act, however, in the full space. As long as one 
uses soft potentials, such as the Viowk^ SRG, UCOM or JISP, convergent NCSM results 
can be obtained. Such NCSM calculations are variational with the HO frequency and 
the basis truncation parameter Amax acting as variational parameters. 

However, the situation is different when standard NN potentials that generate 
strong short-range correlations, such as AVIS, CD-Bonn 2000, and INOY, are used, 
or when a not-large-enough Amax truncation can be reached with the chiral N^LO NN 
potential (in particular, when it is used in combination with the chiral NNN interaction). 
In order to obtain meaningful results in the truncated (or model) space, the inter-nucleon 
interactions need to be renormalized. We need to construct an effective Hamiltonian 
with the inter-nucleon interactions replaced by effective interactions. By meaningful 
results we understand results as close as possible to the full space exact results for a 
subset of eigenstates. Mathematically we can construct an effective Hamiltonian that 
exactly reproduces the full space results for a subset of eigenstates. In practice, we 
cannot in general construct this exact effective Hamiltonian for the A-nucleon problem 
we want to solve. However, we can construct an effective Hamiltonian that is exact for 
a two-nucleon system or for a three-nucleon system or even for a four-nucleon system. 
The corresponding effective interactions can then be used in the 74-nucleon calculations. 
Their use, in general, improves the convergence of the problem to the exact full space 
result with the increase of the basis size. By construction, these effective interactions 
converge to the full-space inter-nucleon interactions, therefore, guaranteeing convergence 
to the exact solution, when the basis size approaches the infinite full space. 

In our approach we employ the so-called Okubo or Lee-Suzuki similarity 
transformation method [231 1211 HO], which yields a starting-energy independent 
hermitian effective interaction. We first recapitulate general formulation and basic 
results of this method. Applications of this method for computation of two- or three- 
body effective interactions are described afterwards. 



2.3.1. Lee-Suzuki similarity transformation method Let us consider an arbitrary 
Hamiltonian H with the eigensystem E^, \ k), i.e., 

H\k) = Ek\k) . (14) 
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Let us further divide the full space into the model space defined by a projector P and the 
complementary space defined by a projector Q, P+Q = 1. A similarity transformation of 
the Hamiltonian e~'^He'^ can be introduced with a transformation operator u satisfying 
the condition u = QuoP. The transformation operator is then determined from the 
requirement of decoupling of the Q-space and the model space as follows 

Qe-'^He'^P = . (15) 

Using a Feshbach construction, one can show that the particular choice of the decoupling 
condition (flSl) ensures that the effective Hamiltonian is energy independent [12]. If we 
denote the model space basis states as \ap), and those which belong to the Q-space, 
as lag), then the relation Qe"^ He'^ P\k) = 0, following from Eq. ( fTSl) . will be satisfied 
for a particular eigenvector \k) of the Hamiltonian (JT^ . if its Q-space components 
can be expressed as a combination of its P-space components with the help of the 
transformation operator u, i.e., 

{aqlk) = '^{aQ\uj\ap){ap\k) . (16) 

ap 

If the dimension of the model space is dp, we may choose a set /C of dp eigenevectors, 
for which the relation f|T6l) will be satisfied. Under the condition that the dpxdp matrix 
defined by the matrix elements {ap\k) for G /C is invertible, the operator uj can be 
determined from ( fT6|) as 

{aQ\uj\ap) = J2Mk){k\ap) , (17) 

where we denote by tilde the inverted matrix of {ap\k), e.g., ^q,^ (fc|«p) («p|fc') = Sk^k', 
for k, k' G /C. 

The hermitian effective Hamiltonian defined on the model space P is then given by 

m 

H,s = [P{1 + uj^uj)P] PH{P + QuoP) [P{1 + Juj)P] . (18) 

By making use of the properties of the operator uo, the effective Hamiltonian can 
be rewritten in an explicitly hermitian form as 

i^eff = [P{1 + uj^uj)P] (P + PJQ)H{QujP + P) 

X [P{1 + uo^uo)PY^'^ . (19) 

With the help of the solution for uj (fT71) we obtain a simple expression for the matrix 
elements of the effective Hamiltonian 

{ap\H^^\ap>) = ^ {ap\k"){}z"\ ap" ) 

Opi, Op, II kk'k"£JC 

X {apn\k)Ek{k\ap'"){ap'"\k'){k'\ap') . (20) 

with all the summations over the Q-space basis states removed. The effective 
Hamiltonian (!20l) reproduces the eigenenergies Ek, k E K. in the model space. 
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It has been shown [?T] that the hermitian effective Hamiltonian fll9p can be obtained 
directly by a unitary transformation of the original Hamiltonian: 

^eff = Pe-'He'P , (21) 

with an anti-hermitian operator S = arctanh(co' — cj^). The transformed Hamiltonian 
then satisfies decoupling conditions Qe~^ He^ P = Pe'^He^Q = 0. 

We can see from Eqs. (l20l) that in order to construct the effective Hamiltonian we 
need to know a subset of exact eigenvalues and model space projections of a subset of 
exact eigenvectors. This may suggest that the method is rather impractical. Also, it 
follows from Eq. (120|) that the effective Hamiltonian contains many-body terms, in fact 
for an A-nucleon system, all terms up to A-hodj will in general appear in the effective 
Hamiltonian, even if the original Hamiltonian consisted of just two-body or two- plus 
three-body terms. 



2.3.2. Two-body effective interaction in the NCSM In the ah initio NCSM we use the 
above effective interaction theory as follows. Since the two-body part dominates the A- 
nucleon Hamiltonian ([2]), it is reasonable to expect that a two-body effective interaction 
that takes into account full space two-nucleon correlations would be the most important 
part of the exact effective interaction. If the NNN interaction is taken into account, 
a three-body effective interaction that takes into account full space three-nucleon 
correlations would be a good approximation to the exact A-body effective interaction. 
We construct the two-body or three-body effective interaction by application of the 
above described Lee-Suzuki procedure to a two-nucleon or three-nucleon system. The 
resulting effective interaction is then exact for the two- or three-nucleon system. It is 
an approximation of the exact A-nucleon effective interaction. 

Using the notation of Eq.([2]), the two-nucleon effective interaction is obtained as 

^2efT,12 = P2[e~'''{hl + h2 + \/i^'^)e^- - {h, + h2)]P2 , (22) 

with 5*12 = arctanh(co'i2 — ci^lg) and P2 is a two-nucleon model space projector. The 
two-nucleon model space is defined by a truncation A^i2max corresponding to the A- 
nucleon A^max- For example, for A = 3,4, A'i2max = A'max, for p-shell nuclei with A> b 
A^i2max = Ainax + 2. The Operator uji2 is obtained with the help of Eq. (fT7|) from exact 
solutions of the Hamiltonian hi -\- h2 + V^'^, which are straightforward to find. In 
practice, we actually do not need to calculate uu, rather we apply Eqs. ( l20l) with the 
two-nucleon solutions to directly calculate P2e~^'^'^ {hi -\- h2 + V^''^)e^^^ P2. To be explicit, 
the two-nucleon calculation is done with 

= Ho2 + Vl",^ = ^ + Imn'r^ + VMV2r) - ^r"^ , (23) 

where r = \J\ifi — ^2) and p = \J\{'Pi — P2) and where Hq2 differs from hi + h2 by 

the omission of the center-of-mass HO term of nucleons 1 and 2. Since V^'^ acts on 
relative coordinate, the Su is independent of the two-nucleon center of mass and the 
two-nucleon center-of-mass Hamiltonian cancels out in Eq. (122|) . We can see that for 
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A > 2 the solutions of (125]) are bound. The relative-coordinate two-nucleon HO states 
used in the calculation are characterized by quantum numbers \nlsjt) with the radial 
and orbital HO quantum numbers corresponding to coordinate r and momentum p. 
Typically, we solve the two-nucleon Hamiltonian (123|) for all two-nucleon channels up 
to j = 8. For the channels with higher j only the kinetic-energy term is used in the 
many-nucleon calculation. The model space P2 is defined by the maximal number of 
allowed HO excitations A''i2max from the condition 2n + l < A''i2max- In order to construct 
the operator uj ( TTTl) we need to select the set of eigenvectors /C. We select the lowest 
states obtained in each channel. It turns out that these states also have the largest 
overlap with the P2 model space. Their number is given by the number of basis states 
satisfying 2n + l < iVi2max- 

The two-body effective Hamiltonian used in the A-nucleon calculation is then 

A A 

HleS = E + E ^^eff., ■ (24) 
i=l i<j 

At this point we also subtract the Hcm and, if the Slater determinant basis is to be used, 
we add the Lawson projection term P{Hcm — ^hfl) to shift the spurious CM excitations. 
Eigenenergies of physical states are independent of the parameter j3. 

2.3.3. Three-body effective interaction in the NCSM An improvement over the two- 
body effective interaction approximation is the use of the three-body effective interaction 
that takes into account the full space three-nucleon correlations. If the NNN interaction 
is included, the three-body effective interaction approximation is rather essential for 
A > 3 systems. First, let us consider the case with no NNN interaction. The three- 
body effective interaction can be calculated as 



T/ININ _ p 

''3efr,123 ~ 3 



-[hi + h2 + h^)]Pz. (25) 

Here, = arctanh(co'i23 — ci;|23) and P3 is a three-nucleon model space projector. 
The P3 space contains all three-nucleon states up to the highest possible three-nucleon 
excitation, which can be found in the P space of the A-nucleon system. For example, for 
A = 6 and A/'max = 6 {6hfl) space we have P3 defined by A'i23max = 8. Similarly, for the 
p-shell nuclei with A> 7 and A^max = 6 {6hQ) space we have A^i23max = 9. The operator 
Ci;i23 is obtained with the help of Eq. (fTTI) from exact solutions of the Hamiltonian 
hi + h2 + + V^'^ + V^'^ -\- V^'^, which are found using the antisymmetrized three- 
nucleon Jacobi coordinate HO basis. In practice, we again do not need to calculate a;i23, 
rather we apply Eqs. (l20l) with the three-nucleon solutions. The three-body effective 
interaction is then used in A-nucleon calculations using the effective Hamiltonian 

Ht. = J:h^ + J^ E ^3X..' (26) 

i=l i<j<k 

where the factor takes care of over-counting the contribution from the two-nucleon 
interaction. 
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If the NNN interaction is included, we need to calculate in addition to the 
following effective interaction 



T/NN+NNN _ p 



e "^123 {^hi + h2 + hs + V^2 + ^^13 + ^23' 



+ H^NN,123)e^i23 _ ^hi + h2 + /l3)J ^3 • (27) 

This three-body effective interaction is obtained using full space solutions of the 
Hamiltonian hi -\- h2 + -\- V^'^ + V^'"^ + ^23'"^ + Vnnn,i23- The three-body effective 
interaction contribution from the NNN interaction we then define as 

T/NNN _ T/NN+NNN _ t/NN fr.Q\ 
^3eff,123 — '^3efr,123 '^3ofr,123 ■ K"^") 

The three-body effective Hamiltonian used in the A-nucleon calculation is then 

1=1 i<j<k i<j<k 

As in the case of the two-body effective Hamiltonian fl2^ . we subtract the -ffcM and, 
if the Slater determinant basis is to be used, we add the Lawson projection term 

It should be noted that all the effective interaction calculations are performed in 
the Jacobi coordinate HO basis. As discussed above, the two-body effective interaction 
is performed in the \nlsjt) basis and the three-body effective interaction in the \NiJT) 
basis f|TT]) . In order to perform the A-nucleon calculation in the Slater determinant HO 
basis, as is typically done for A > 4, the effective interaction needs to be transformed 
to the single-nucleon HO basis. This is done with help of the HO wave function 
transformations. The details for the three-body case, in particular, are given in Refs. [13] 
and [44j . 

It should also be noted that one may attempt to separate the two-body and the 
three-body parts of the ( [25l) . This has not been done yet in the NCSM calculations 
as the current implementation ( 126|) proved robust (as also demonstrated in the next 
section). In recent one-dimensional model calculations with SRG evolved interactions 
such a separation has been achieved and shown to be useful [l^. It should be also 
explored within the NCSM, although care must be taken to avoid introducing spurious 
model-space effects. 



2.4- Effective operators 

Besides spectra, other properties of the nuclear states are of interest, as they impose 
a strong test on the theoretical wave functions. For consistency, the same unitary 
transformation used to compute the effective interaction should be employed in order 
to obtain effective operators in the model spaces used to diagonalize the effective 
Hamiltonian. 

In addition to consistency, another motivation for implementing the renormalization 
of general operators is the long standing effective charge problem in the 
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phenomenological shell model. Arising from the inevitable truncation of the Hilbert 
space, the relatively large effective charges were found to be essential in the overall 
description of the transition strength. However, previous perturbation theory attempts 
to describe phenomenological charges needed to obtain correct transition strengths have 
been unsuccessful |16]. On the other hand, investigations within the framework of the 
NCSM have reported some progress in explaining the large values of the effective charges 

m- 

The renormalization of effective operators is much more involved than the 
renormalization of the Hamiltonian. In order to ensure energy independence of the 
effective operator, the decoupling condition (JTSll has to be supplemented with the 
Hermitian conjugate |12]. This transformation, however, has the advantage that the 
effective Hamiltonian produced is Hermitian, as discussed in Sec. 12.31 and written out 
explicitly in Eq. ( |T9l) . 

A general tensor operator can change the spin and isospin. Hence, the 
renormalization of a rank A J, AT tensor operator writes as 

q(AJ;AT) _ Pj'T' + Pj'T"^\iT'Qj'T' q{AJ;AT) ^JT + QjT^JtPjT ^^q^ 



PjiT' + u\,j,,UjiT' y PjT + i^jT^JT 



where we have shown explicitly the possible change in spin and isospin of the initial 
(J, T) and final (J', T') states, respectively. Equation fl30l) is the generalization of 
f fT9|) . and shows the complexity of the renormalization of a a tensor operator compared 
with the renormalization of a scalar operator, as the tensor operator allows for the the 
possible mixture of different spin and isospin quantum numbers. Finally, because the 
transformation is a scalar, the effective operator preserves the tensor character of the 
starting operator. 

The simplest approximation for the unitary transformation is at the two-body 
cluster level. Because of the complexity of the renormalization, the two-body cluster 
is the only one developed so far for general operators [HI HHl [50]. Under this 
approximation, the transformation becomes 

A 

S^^Y. (31) 

i>j=i 

with Sij = arctanh(co',j — ujjj). Applying the operator identity 

e'^'Oel = 0+[0, S2] + ^[[O, S^], S^] + ... (32) 
to transform a general one-body operator O^-^^ = J2f=i Oi, one obtains 

A A 

= 0(1) + ^ + o„ S,.j\ + Y,[[0^ + 0„ S,j\,S,j\ + (33) 

i>j=l i>j 

where we have retained only the one- and two-body terms, neglecting higher body 
contributions, such as [Oj, Sjk], with i ^ j and i ^ k. Resummation of the commutators 
yields 
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i 

A 

+ + ^^•) - + ^^■)] ^2- (34) 

i>j=i 

Analogously, for a general two-body operator 

A 

P202esP2 = P2Y1 e'^- 0,,e^- P2 , (35) 
i>j=l 

and, in particular, the effective Hamiltonian is given by 

A 

P2H2esP2 = P2j2^^P^ 

i=l 

A 

+ [e"^'' + + ^^i) e^'' - - ^i] ^2, (36) 

i>j=i 

which recovers expression fl22l) for the effective interaction. 

^.5. Convergence tests 

In this subsection, we give examples of convergence of ab initio NCSM calculations. 

First, we discuss calculations for s-shell nuclei. In Fig. [H we show the convergence 
of the and ^He ground-state energies with the size of the basis. Thin lines correspond 
to results obtained with the NN interaction only. Thick lines correspond to calculations 
that also include the NNN interaction. Here, we use the chiral effective field theory 
(EFT) NN interaction of Ref. [30] and the local chiral NNN interaction that will 
be discussed in detail in the next section. The solid lines correspond to '^H (^He) 
calculations with two-body (three-body) effective interaction derived from the chiral 
EFT potentials. The dashed lines correspond to calculations with the bare, that is the 
original, unrenormalized chiral EFT interactions. In calculations, the bare NNN 
interaction is added to either the bare NN (dashed thick line) or to the effective NN 
interaction (solid thick line). We observe that the convergence is faster when the effective 
interactions are used. However, starting at about A^max = 24(18) the convergence is 
reached in ^H(^He) calculations also with the bare NN interaction. It should be noted, 
however, that p-shell calculations with the NNN interactions are presently feasible in 
model spaces up to A'max = 6 or A'max = 8. The use of the three-body effective interaction 
is then essential in the p-shell calculations with NN+NNN interactions. It should be 
noted that in calculations with the effective interaction, the effective Hamiltonian is 
different at each point, as the effective interaction depends on the size of the model space 
given by A^max- The calculation with the bare interaction is a variational calculation 
converging from above with Amax and HO frequency fl as variational parameters. The 
calculation with the effective interaction is not variational. The convergence can be 
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Figure 1. '^H (left) and ''He (right) ground-state energy dependence on tiie size of 
the basis. The HO frequencies of HQ = 28 MeV (^H) and 28 or 36 MeV (''He) were 
employed. Results with (thick lines) and without (thin lines) the NNN interaction are 
shown. The solid lines correspond to calculations with two-body (^H) or three-body 
^He) effective interactions, the dashed lines to calculations with the bare interactions. 



from above, from below or oscillatory. This is because a part of the exact effective 
Hamiltonian is omitted. The calculation without NNN interaction converges to the 
ground-state energy —7.852(5) MeV, well above the experimental —8.482 MeV. Once the 
NNN interaction is added, we obtain —8.473(5) MeV, close to experiment. As discussed 
in the next section, the NNN parameters were tuned to reproduce the average of the 
and ^He binding energies. 

The rate of convergence also depends on the choice of the HO frequency. The "^He 
calculations without NNN interaction were done for two different HO frequencies. It 
is apparent that convergence to the same result occurs in both cases. We note that 
in the case of no NNN interaction, we may use just the two-body effective interaction 
(two-body cluster approximation), which is much simpler. The convergence is slower, 
however, see discussion in Ref. [53|- We also note that ^He properties with the chiral 
EFT NN interaction that we employ here were calculated using the two-body cluster 
approximation in Ref. [55] and present results are in agreement with results found 
there. Our ^He ground-state energy results are —25.39(1) MeV in the NN case and 
-28.34(2) MeV in the NN-^NNN case. The experimental value is -28.296 MeV. 
We note that the present ab initio NCSM ^H and ^He results obtained with the 
chiral EFT NN interaction are in a perfect agreement with results obtained using 
the variational calculations in the hyperspherical harmonics basis as well as with the 
Faddeev-Yakubovsky calculations published in Ref. ^6]. A satisfying feature of the 
present NCSM calculation is the fact that the rate of convergence is not affected in any 
significant way by inclusion of the NNN interaction. 

As an example of convergence of ab initio NCSM calculations for p-shell nuclei, we 
present ^Li results obtained using the INOY and the chiral EFT NN potential. The 
dependence of the NCSM absolute and excitation energies on the basis size is presented 
in Fig. [2j The calculations were performed using the two-body effective interaction in 
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the Slater determinant HO basis with the shell- model code Antoine [STj. Results for 
other HO frequencies were published in Refs. |55l [58]. As discussed in Ref. [55], the 
convergence rate with iVmax is different for different states. In particular, the 3"'"0 state 
and the 0"'"1 state converge faster in the higher frequency calculations {HQ = 12, 13 MeV 
with the chiral EFT NN potential), while the higher lying states converge faster in the 
lower frequency calculations [hfl = 8, 10 MeV with the chiral EFT NN potential). The 
results on the right of Fig. [2] demonstrate a good convergence of the excitation energies, 
in particular, for the 3"'"0 and 0'''1 states. An interesting result is the overestimation of 
the 3"'"0 excitation energy compared to experiment, in particular with the chiral EFT 
NN potential. It turns out that this problem is resolved once the NNN interaction is 
included in the Hamiltonian. More discussion on eigenenergy convergence in p-shell 
nuclei NCSM calculations can be found, e.g., in Refs. [59l [60t [6T]. 




Figure 2. Calculated absolute (left) and excitation (right) energies of ^Li obtained in 
0hn-12hil {Ohft-lAhn) basis spaces using two-body effective interactions derived from 
the INOY (chiral EFT) NN potential compared to experiment. The HO frequencies of 
hfl = 14 MeV (left) and 12 MeV (right) were used. 

In the final part of this section, we discuss the renormalization of different 
observables, such as electro-magnetic moments, radii, etc. The first exclusive 
investigation of the effective operators in the NCSM framework concentrated on 
implementing a procedure in which the relative states in the excluded space were 
restricted, and observed convergence by releasing the restriction [18] . The goal of the 
investigation in Ref. [H] was a test of the numerical implementation in a simple model, 
in which the "full" space was restricted to a numerically tractable size, so that all 
missing correlations could be exactly calculated. Even in this simple model it was found 
that the two-body cluster renormalization was very weak for long range operators, such 
as E2 transitions. The same results were later reported when the same method was 
implemented in realistic cases |49j . 

The main goal of Ref. |19] was a qualitative understanding of the influence 
of effective operators and not a highly accurate description of the experimental 
data; therefore, the NNN interactions were left out. The same goal also motivated 
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the use of rather small model spaces and of the two-body cluster approximation, 
given that, as expected from the convergence properties of effective operators, larger 
renormalization effects are expected in smaller model spaces. Furthermore, a more 
efficient implementation of the renormalization than the one in Ref. [H] was introduced, 
similar to the one applied to the Hamiltonian. Thus, the renormalization was 
implemented in relative coordinates, allowing the same treatment for general operators 
as for the Hamiltonian, as long as the former can be written in relative coordinates |49j . 
This implementation was tested on the deuteron, where the two-body cluster provides 
the exact solution. In that case, the bare quadrupole operator in 4hQ gave 0.179 e 
fm^ for the quadrupole moment, while the value of 0.270 e fm^, described by the AV8' 
potential, was obtained using the corresponding effective operator in the same model 
space. However, when the same procedure has been applied to a realistic many-body 
problem, the result was different. A very weak renormalized B{E2) value was obtained 
even in small model spaces. As an example, we have looked at the B{E2; 3"'"1 — » 1^1) in 
^Li, where a 2hQ calculation gives 2.647 fm^ when the bare operator is employed, and 
2.784 fm^ when the effective operator is used [19]. These results have been obtained 
with the Argonne V8' NN interaction. A calculation for the same observable, but with 
the CD-Bonn 2000 NN interaction, which is expected to give comparable results with 
the AV8' potential, obtained B{E2; 3+1 1+1) = 10.221 fm^ with the bare operator 
in lOhfl model space. Overall, the difference between the bare operator results in the 
2Hfl and lOHQ model spaces, coupled with the small renormalization at the two-body 
cluster level, indicate sizable effective many-body effects needed to correct the 2hQ 
B{E2) value. 

The case of the kinetic energy operator is completely different from the quadrupole 
transitions presented above. In Refs. pTl B9] . large renormalization was obtained even 
at the two-body cluster for the kinetic energy. The kinetic energy is short range, while 
the quadrupole is long range; at the two-body cluster level, the unitary transformation 
renormalizes mainly the short-range core of the interaction, leaving unchanged the long 
range part. Hence, in order to account for long-range correlations in the two-body 
cluster approximation, one needs to enlarge the model space. To test this hypothesis, 
a Gaussian operator of variable range was used in Ref. ^9]. There, by observing the 
variation with the model space / HO frequency of the expectation values calculated with 
the bare and effective operators for several different ranges, it was demonstrated that 
a short-range two-body operator is renormalized accurately at the two-body cluster 
level, while a long-range operator is weakly renormalized. To further illustrate the 
power of the unitary transformation approach to the renormalization of short-range 
operators, we turn to an observable probing short-range correlations. The inclusive 
(e, e') longitudinal data presents one of the clearest experimental signatures for short- 
range correlations in the wave-function of the ground state, at least for light nuclei. 
A quantitative measure of the short-range correlations is the longitudinal-longitudinal 
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Figure 3. The longitudinal-longitudinal distribution function in ''He, obtained using 
bare operators (left panel) and effective operators (right panel). The HO energy used 
in this calculation was hfl — 28 MeV, while the NN interaction was CD Bonn. 

distribution function (connected to the Coulomb sum rule) [52l [53] 
PLL(g) = ^{g-S.\jo{q\ri - rj\){l + r^,i)(l + T,j)\g.s.), 

where Jq is the spherical Bessel function of zero order, and q the momentum transfer. 

In Figure [3] we present the results for the longitudinal-longitudinal distribution 
function for ^He. At high momentum transfer, the results obtained using bare operators 
depend strongly upon the model space. On the other hand, the results obtained with 
effective operators are model space invariant at high q, although in these model spaces 
the wave function is not fully converged, since the energy is not converged in these very 
small model spaces (see Fig. 3 in Ref. [50]). They agree with the values computed in 
larger model spaces. At intermediate momentum transfer, i.e., q ~ 2.5 fm~^, even the 
effective operator results vary. This effect is due to the fact that the long range part of 
the operator has not yet converged in these small model spaces. In larger model spaces, 
where the long-range correlations are better described, the agreement is even better. 

Similar results for the longitudinal-longitudinal distribution function have been 
obtained for ^^C, where calculations in very large model spaces are not possible. 
However, even in the smallest model space, Ohfl, we were able to obtain good results 
for high momentum transfer, which reproduce the values in larger model spaces j50j . 

In conclusion, short-range operators (high momentum transfer) are very well 
renormalized and the results become model-space independent even in the two- 
body cluster approximation, while long-range operators, such as the quadrupole 
transition operator, or the longitudinal-longitudinal distribution function for small and 
intermediate momentum transfer, are only weakly renormalized. It all comes down to 
the effects of the unitary transformation, which, as discussed before, at the two-body 
level renormalizes only the short-range part of the interaction, while the long-range part 
is recovered in larger model spaces. 
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3. Light nuclei from chiral EFT interactions 

Interactions among nucleons are governed by quantum chromodynamics (QCD). In 
the low-energy regime relevant to nuclear structure, QCD is non-perturbative, and, 
therefore, hard to solve. Thus, theory has been forced to resort to models for the 
interaction, which have limited physical basis. New theoretical developments, however, 
allow us connect QCD with low-energy nuclear physics. The chiral effective field 
theory (xEFT) [62] provides a promising bridge. Beginning with the pionic or the 
nucleon-pion system [63] one works consistently with systems of increasing nucleon 
number [HH [65l [66]. One makes use of spontaneous breaking of chiral symmetry to 
systematically expand the strong interaction in terms of a generic small momentum 
and takes the explicit breaking of chiral symmetry into account by expanding in the 
pion mass. Thereby, the NN interaction, the NNN interaction and also vrN scattering 
are related to each other. The xEFT predicts, along with the NN interaction at the 
leading order, an NNN interaction at the third order (next-to-next-to-leading order or 
N^LO) [621 [671 [68], and even an NNNN interaction at the fourth order (N^LO) [69]. The 
details of QCD dynamics are contained in parameters, low-energy constants (LECs), not 
fixed by the symmetry. These parameters can be constrained by experiment. At present, 
high-quality NN potentials have been determined at order N^LO [30]. A crucial feature 
of xEFT is the consistency between the NN, NNN and NNNN parts. This consistency 
also extends to the nuclear current. As a consequence, at N^LO and N^LO, except for 
two LECs, assigned to two NNN diagrams, the potential is fully constrained by the 
parameters defining the NN interaction. 

We adopt the potentials of the xEFT at the orders presently available, the NN at 
N^LO of Ref. [30] and the NNN interaction at N^LO [671 [68]. Since the NN interaction 
is non-local, the ab initio NCSM is the only approach currently available to solve the 
resulting many-body Schrodinger equation for mid-p-shell nuclei. We are in a position 
to use the ab initio NCSM calculations in two ways. One of them is the determination 
of the LECs assigned to two NNN diagrams that must be determined in A > 3 systems. 
The other is testing predictions of the chiral NN and NNN interactions for light nuclei. 

3.1. Chiral ISPLO three-nucleon interaction 

The NNN interaction at N^LO of the xEFT is comprised of three parts: (i) The two-pion 
exchange, (ii) the one-pion exchange plus contact and (iii) the three-nucleon contact, 
see Fig. HI In this work, we regulate the the NNN terms with a regulator depending 
on the momentum transfer similarly as done, e.g., for the Tucson-Melbourne NNN 
interaction [70j, which results in a local xEFT NNN interaction. This is advantageous 
for some many-body approaches, including the NCSM, because a local NNN interaction, 
in particular the two-pion exchange term, is easier to implement. Full technical details 
are given in Ref. [7T]. 

The LECs associated with the two-pion exchange also appear in the NN interaction 
and are, therefore, determined in the A = 2 system. The one-pion exchange plus contact 
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term (D-term) is associated with the LEG and the three-nucleon contact term (E- 
term) is associated with the LEG ce- It is interesting to note that cd represents 
a renormahzation of the short range interaction of two nucleons, rather than three. 
Gonsequently, it manifests itself not only in the A^A^ — it — N contact term of the 
interaction, but also in the two-nucleon contact vertex with an external probe of the 
exchange currents. 

The Cd and ce LEGs, expected to be of order one, can be constrained by the A = 3 
binding energy. One then still needs an additional observable to determine the two 
parameters. The first determination of c^i and ce was attempted using as constraints 
the binding energy and nd doublet scattering length, and adopting the full interaction 
up to N^LO [68j. However, this proved to be difficult due to a correlation between these 
two observables, and the large experimental uncertainty on the scattering length. Later, 
the N^LO NN potential was combined with the NNN at N^LO (non-local, regulated with 
nucleon momenta) to study the "^Li structure [S]. In this work, besides the binding 
energy the second constraint on the undetermined LEGs was the energy of the ^He 
ground state. As a result of the correlation between these two observables, known as 
Tjon line, fitting the ground-state energy automatically results in a ^He binding 
energy that is within a few hundred keV of experiment. The subsequent fine-tuning of 
this binding energy is then very sensitive to the structure of the adopted NNN force. 
Hence, small variations of the cutoff, different regularization schemes, missing terms of 
the interaction, etc., tend to produce large swings in the extracted values of cd and ce- 
A different approach that we describe and expand on it here was adopted in Ref. |72j . 
There, a preferred choice for the two LEGs was obtained by complementing the constrain 
on the A = 3 binding energies with a sensitivity study on the radius of the '^He and on 
various properties of p- shell nuclei. 

Fig. [5] shows the trajectories of the two LEGs cd and ce that were determined 
in Ref. [72] from fitting the binding energies of the A = 3 systems. Separate curves 
are shown for ^H and ^He fits, as well as their average. We also show the calculated 
ground-state energy and charge radius of ^He obtained with the constrained LEGs. As 
demonstrated in Fig. [1] in Section 12.51 our A = 3 and A = 4 results presented in Fig. [5] 
are fully converged. There are two points where the binding of ^He is reproduced exactly. 
One of them for c/j ~ 1 and the other with c/j ~ 10. However, as a consequence of 
the correlation bewteen the triton and '^He binding energies, we observe that in the 
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Figure 5. Relations between cb and ce for which the binding energy of '^H (8.482 
MeV) and "^He (7.718 MeV) are reproduced, (a) "^He ground-state energy along the 
averaged curve, (b) "^He charge radius Tc along the averaged curve. Dotted lines 
represent the uncertainty due to the uncertainties in the proton charge radius. 

whole investigated range oi cr, and c^, the calculated ^He binding energy is within a 
few hundred keV of experiment. Consequently, the determination of the LECs in this 
way is likely not very stringent. By examining calculations of the ^He charge radius, we 
can see that a result consistent with experiment, taking into account the uncertainty of 
the proton charge radius, is obtained for c^) values in the range from ~ — 2 to ~ +2. 
This observation rules out the range of the large cb ~ 10 values that overestimate the 
^He radius and on top of it might be considered "unnatural" from the xEFT point of 
view. 

In Tables [H and [21 we present a collection of A = 3 and 4 data, respectively, 
obtained with and without inclusion of the NNN force for = 1 [ce = —0.029). 
The corresponding ground-state energy convergence within the NCSM was shown in 
Figs. [1] in Section 12.51 Besides the triton ground-state energy, which is by construction 
within a few keV of experiment, the NN-I-NNN results for the nd doublet and quartet 
scattering lengths and ''He ground-state energy and point-proton radius are in perfect 
agreement with measurement. While for the n^H singlet scattering length the inclusion 
of the NNN force worsens the disagreement with respect to experiment to some extent, 
the n^H triplet scattering length improves with the NNN included. We also note a 
perfect agreement between the two theoretical approaches, the ab initio NCSM and the 
variational HH method of Ref . [I^ , for the bound-state results. 

As the €£) LEC enters also the nuclear current, it is possible to utilize, e.g., the 
triton half life, as another observable in addition to the A = 3 binding energy constraint 
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Table 1. "^H g.s. energies (in MeV), point-proton radii (in fni) and nd scattering 
lengths (in fm), obtained using the N^LO NN potential 30j with and without the local 
N^LO NNN interaction ^Tlj with cd = 1 and ce = —0.029, compared to experiment. 
Calculations performed within the NCSM and/or hyperspherical harmonics (HH) 
expansion approaches. 

nd 



NN NCSM m] -7.852(5) 1.650(5) 

NN HH [n] -7.854 1.655 1.100 6.342 

NN+NNN NCSM [H] -8.473(5) 1.608(5) 

NN+NNN HH [17] -8.474 1.611 0.675 6.342 

Expt. -8.482 1.60 

Expt. (Tuna [76] - - 0.65(4) 6.35(2) 

Expt. [771 [78] - - 0.645(8) 



Table 2. The same as in Table [T]for ''He and n^H. 
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for the determination of the NNN LECs. This was done recently in Ref. [73] . Taking 
into account experimental errors, a narrow range of cd values around co ~ —0.2 (the 
corresponding ce = —0.205 from the A = 3 binding energy) was found to agree with 
the measurements. An interesting observation was made in Ref. [T3] , namely, the triton 
half life is quite insensitive, unlike most other observables, to the NNN terms in the 
Hamiltonian. This is makes the LEC determination of Ref. [T3] rather robust. 

3.2. Results for p- shell nuclei 

Sensitivity of the p-shell nuclear properties to the choice of the cd and ce LECs was 
investigated in Ref. [T2]. First, the A = 3 binding energy constraint was maintained. 
Second, the sensitivity study was limited to the cd values in the vicinity of the point 
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Figure 6. Dependence of the ^Li quadrupole moment (left) and the B(E2;3]^0 — > 
lj^0)/B(E2;3|j''0 ^2^) ratio (right) on the value of cjj with the constrained by 
the A = 3 binding energy fit for different basis sizes. The HO frequencies of hil = 13 
MeV (left) and 14 MeV (right) were employed. 



cd ~ 1, in particular values from —3 to +2 that include the range compatible with the 
^He radius. 

While most of the p-shell nuclear properties, e.g., excitation spectra, are not very 
sensitive to variations of c/j in the vicinity of the c^i ~ 1 point, we were able to identify 
several observables that do demonstrate strong dependence on cn. For example, the 
^Li quadrupole moment changes sign depending on the choice of co, as can be seen in 
Fig. El In the right of Fig. [6], we display the ratio of the B(E2) transitions from the 
ground state to the first and the second 1^0 state. This ratio changes by several orders 
of magnitude depending on the cd variation. This is due to the fact that the structure 
of the two 1"''0 states is exchanged depending on c^i. In addition, cd dependence of the 
^^C B(M1) transition from the ground state to the 1+1 state was discussed in Ref. [72] . 
Also, the importance of the NNN interaction in reproducing the experimental value was 
illustrated [721 [79]. Overall the results show that for cd < —2 the '^He radius and the 
^Li quadrupole moment underestimate experiment, while for cd > the lowest two 1+ 
states of are reversed and the ^^C B(M1;0"*'0 —>■ 1+1) is overestimated. Therefore, 
the value of cd = — 1 was chosen in Ref. [72] as globally the best choice. We note that 
the triton half-life study suggests a cc value in the range around cd ~ —0.2. This 
result is not inconsistent with the study of the p-shell nuclei. It is straightforward to 
reconcile these findings. First, one may consider the p-shell calculations less reliable 
than the much-less-involved A = 3 calculations. However, it is quite plausible that 
the re-normalization of the cd value for p-shell nuclei mimics the effect of (neglected) 
higher-order NNN force terms, which are irrelevant for the calculation of the triton 
half life. In fact, a closer look at the ^He results shown in Fig. [5l and Table [21 the ^Li 
results from Fig. [6l the ^''B results from the right of Fig. [6] and the ^^C B(M1) results, 
suggests a drift of the optimal value towards smaller (increasingly more negative) 
values with nuclear mass. It is natural to expect that an effect of the higher order NNN 
terms will become more important with increasing mass of the nucleus. Another issue 
that deserves attention is the determination of the C3 and C4 LECs from the NN data 
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Figure 7. -'^^B excitation spectra as function of the basis-space size iVmax with the 
chiral NN (left) and NN+NNN (right) interaction compared to experiment. The isospin 
of the states not exphcitly depicted is T=0. The HO frequencies of hil — 14 MeV (left) 
and 15 MeV (right) were used. 



and the assessment of the extent to which they would influence the determination of cd 
from the triton half life and the p-shell nuclei calculations. The C4 LEG is, in particular, 
poorly constrained by the NN data fit pO]. 

In Figs. O we present the excitation spectra of ^°B, as a function of A^max, for both 
the chiral NN+NNN, as well as with the chiral NN interaction alone. In both cases, the 
convergence with increasing A^max is quite reasonable for the low-lying states. Similar 
convergence rates are obtained for our other p— shell nuclei calculations. 

A remarkable feature of the results is the observation that the chiral NN 
interaction alone predicts the incorrect ground-state spin of ^'^B. The experimental value 
is 3"'"0, while the calculated one is 1"'"0. On the other hand, once we also include the 
chiral NNN interaction in the Hamiltonian, which is actually required by the xEFT, the 
correct ground-state spin is predicted. Further, once we select the cd value, as discussed 
above, i.e., cd = —1, we also obtain the two lowest 1"'"0 states in the experimental order. 

We display in Fig. [8] the natural-parity excitation spectra of four nuclei in the 
middle of the p— shell with both the NN and the NN+NNN effective interactions from 
xEFT. The results shown are obtained in the largest basis spaces achieved to date for 
these nuclei with the NNN interactions, A^max = 6 {6hQ). Overall, the NNN interaction 
contributes significantly to improve theory in comparison with experiment. This is 
especially well-demonstrated in the odd mass nuclei for the lowest, few excited states. 
The case of the ground-state spin of ^°B and its sensitivity to the presence of the NNN 
interaction, discussed also in Fig. [71 is clearly evident. We note that the ^°B results in the 
left panel of Fig. U\ only with the NN interaction, were obtained with the HO frequency 
of HQ = 15 MeV, while those in Fig. [8] are with HQ = 14 MeV. A weak HO frequency 
dependence of the iV^ax = 6 results is evident. The ^*^B results with the NN+NNN 
interaction, presented in Figs. [7] and [HI were obtained using the same HO frequency. 
Still, one may notice small differences of the A'max = 6 results. The reason behind those 
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Figure 8. States dominated by p-sliell configurations for ^"B, ^^B, ^^C, and ^^C 
calculated at A'^max = 6 using hft = 15 MeV (14 MeV for ^°B). Most of the eigenstates 
are isospin r=0 or 1/2, the isospin label is explicitly sho"wn only for states "with T—1 
or 3/2. The excitation energy scales are in MeV. 



differences is the use of two alternative D-term regularizations, discussed in Ref . [7T] . As 
the dependence on the regulator is a higher-order effect than the xEFT expansion order 
used to derive the NNN interaction, these differences should have only minor effect. It 
is satisfying that the present ^''B results appear to support this expectation. 

Concerning the ^^C results, there is an initial indication that the chiral NNN 
interaction is some'what over-correcting the inadequacies of the NN interaction since, 
e.g., the 1+0 and 4+0 states in ^^C are not only interchanged, but they are also spread 
apart more than the experimentally observed separation. In the ^^C results, we can 
also identify an indication of an overly strong correction arising from the chiral NNN 
interaction, as seen in the upward shift of the | state. However, the experimental | 
may have significant intruder components and is not well-matched with our state. In 
addition, convergence for some higher-lying states is affected by incomplete treatment 
of clustering in the NCSM. This point will be elaborated upon later. These results 
required substantial computer resources. A typical A^max = 6 spectrum, shown in Fig. 
[HI and a set of additional experimental observables took 4 hours on 3500 processors of 
the LLNL's Thunder machine. The A-nucleon calculations were performed in the Slater 
determinant HO basis using the shell model code MFD [81] . 

The calculations presented in this section demonstrate that the chiral NNN 
interaction makes substantial contributions to improving the spectra and other 
observables. However, there is room for further improvement in comparison with 
experiment. In these calculations we used a strength of the 27r-exchange piece of the 
NNN interaction, which is consistent with the NN interaction that we employed (i.e. 
from Ref. [30]). As we already discussed, this strength is somewhat uncertain (see, e.g., 
Ref. [S]). Therefore, it will be important to study the sensitivity of our results with 
respect to this strength. Further on, it will be interesting to incorporate sub-leading 
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NNN interaction terms that are currently under development [H2] and also four-nucleon 
interactions, which are also order N^LO [69]. Finally, it is useful and currently feasible to 
extend the basis spaces to A^max = 8 {8hQ) also for A > 6 to further improve convergence. 

3.3. Photodisintegation of^He 

Over the years, measurements of the a-particle photodisintegration in the near-threshold 
region have been controversial, particularly concerning the height of the cross section 
at the peak, for which one can find differences up to a factor of two between different 
experiments (see, e.g., Ref. [83] and references therein). On the other hand, theoretical 
calculations of the ^He photoabsorption cross section have reached an unprecedented 
level of accuracy in the recent past [831 El] • Predictions obtained using high-precision 
NN and NNN interactions models lie in a rather contained band, which is remarkable 
compared to the large discrepancies still present among the different experimental data. 
In this section, we review the results obtained for the ^He total photoabsorption cross 
section using the above presented chiral NN-I-NNN interactions (for co = —1). 

Nuclear photoabsorption processes can be described in good approximation by the 
cross section 



where u is the incident photon energy and R{lj) is the inclusive response function in 
the long wavelength approximation 



This is the sum of all the transitions from the ground state |\E'o) (of energy Eq) to 
the various allowed final states l^^f) (of energy Ef) induced by the electric dipole 
operator D. Here, a fully ab initio result for the response function (hence for the cross 
section) was obtained by means of the Lorentz integral transform (LIT) method [85] 
implemented in the framework of the NCSM approach [86] • More specifically, we first 
performed an accurate NCSM calculation of the ^He ground state reaching the same 
level of convergence as shown in Fig. [1], by means of effective interactions at the three- 
body cluster level. We then evaluated the LIT of the response f l38l) by applying the 
Lanczos algorithm to the chiral Hamiltonian, using as starting vector the transition 
operator acting on the ground state, D\ip()). Indeed the LIT can be expressed in terms 
of a continued fraction of the so-called Lanczos coefficients (i.e., the elements of the 
tridiagonal Hamiltonian in the Lanczos basis) [87]. After inversion of the LIT |88j . 
the cross section is obtained from Eq. (137]) . The use of three-body effective interaction 
allows a stable and accurate convergence not only of the ^He ground state, but also of 
the photoabsorption cross section. This is presented in the left panel of Fig. [9], where 
the inset shows the tail of the cross section. Due to the selection rules induced by 
the dipole operator, the Lanczos vectors have parity opposite to the ^He ground state. 
Therefore the convergence of our results is studied as a function of A^max/A^max + 1, the 
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Figure 9. Tlie '*He photo-absorption cross section as a function of the excitation 
energy ld. Present resuhs are for xEFT interactions, and in particular: (left panel) 
convergence pattern of the NN+NNN calculation with respect to the model-space 
truncation A/'max for Ml — 28 MeV; (right panel) comparison to the most recent 
experiments [89 l [90l [9T] . 



truncations adopted for the 0"'"0/l~l model spaces, respectively. Our results, presented 
in the right panel of Fig. [9] together with the most recent experiments, show a peak 
around the excitation energy oi uj = 27.8 MeV, with a peak height mildly sensitive to 
the NNN force. The experimental situation in the near-threshold region is controversial: 
two direct measurements performed using quasi-mono-energetic photons [89l |90] show 
discrepancies up to a factor of two on the absolute height of the cross-section peak. We 
find an overall good agreement with the photo-disintegration data from bremsstrahlung 
photons of Nilsson et al. [90], while we reach only the last of the experimental points of 
Ref. |89|. In particular, the confused experimental situation drawn by these two data 
sets does not allow to assess the role of the NNN force effect. Recently Nakayama et 
al. performed an indirect measurements of the a-particle total photo-absorption cross 
section [91] by observing its analog via the '^He(^Li,^Be) reaction at an incident energy of 
455 MeV and at forward scattering angles. Although the uncertainty on this extracted 
absolute cross section is 20% or more, the inclusion of the NNN terms of the interaction 
appear to improve the agreement of the calculated cross section with the latter indirect 
measurement. 

4. Radii, moments and transitions in light nuclei 

4-1. Electric dipole moment of^He 

Massive (and expensive) experimental efforts are directed toward high precision tests of 
the Standard Model. While obviously the experimental confirmation of the Higgs boson 
is of extreme importance, somewhat smaller scale experiments can be set up in search 
for physics beyond the standard model. Presently, several experimental programs are 
pushing the limits on the detection of electric dipole moments (EDMs) of the nucleon, 
nuclei and atoms. Thus, a permanent EDM of a system requires a direct violation of 
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the time-reversal (T) and parity (P) and hence CP violation through the CPT theorem. 
The standard model allows for a very little CP-violation, at levels that are too small to 
be observed currently, and thus any non-zero measurement of EDMs would be a clear 
signal of sources of CP-violation beyond the Standard Model. 

Motivated by the recent proposal of a new scheme for measuring EDMs of light 
nuclei stripped of their electrons (e.g., deuteron, ^He) in magnetic storage rings, the 
first comprehensive ab initio calculation of the EDM of ^He has been reported in Ref. 
[92] . In this investigation, two distinct contributions have been considered: (i) the 
intrinsic EDMs of each nucleon, dp and dn respectively, and (ii) polatization effects 
induced by the CP violation in the NN interaction. 

The three-body problem was solved with high accuracy in relative coordinates, 
using several potential models. Thus, we have obtained results for the CD Bonn 
potential |28j, Argonne VIS interaction [27], as well as the latest generation chiral 
two- [30] and three-body interactions [71]. When three-body forces are included, an 
excellent description of the ground-state properties is achieved, as discussed above. The 
polarization effects, which give the largest contribution to the EDM, were obtained in 
perturbation theory, using the Podolsky method [93], similar to the LIT application to 
exclusive processes [87]. However, in the absence of a EFT derivation of the CP- violating 
interaction, a one-boson exchange model had to be used in this case. Thus, we have 
considered all possible vr, p and uj CP- violating exchanges (for a detailed expression of 
the CP- violating interaction used, see Ref. [92]), the final expression of the EDM being 
expressed as a superposition of CP-violating coupling constants G^, where x stands for 
the meson exchanged, and T for the isospin. The coupling constants are unknown, and 
only limits exist. 

A consistent approach would require that the same transformation used to derive 
the effective interaction be used for any observable calculated with the respective wave 
function. In this case, this means that the dipole-moment operator and the CP-violating 
interaction should be also renormalized, as discussed in Sec. 12. 4[ However, as shown 
in Sec. 12.51 the long-range operators, like the dipole transition operator, are insensitive 
to the renormalization. The long-range correlations are built by increasing the model 
space, in this case by increasing the number of HO shells used to construct the many- 
body basis. Hence, we observe the convergence of the EDM with the model space, and, 
thus, we find that in large-sized model spaces, the results become independent of the 
parameters used (HO frequency and number of HO shells). 

Complete results for the tt, p and u exchanges have been reported in Ref. [92]. In 
Fig. [To], we present the convergence of the isoscalar, isovector and isotensor components 
of the EDM for the vr and p exchanges {u exchange has the same order of magnitude 
as the p). In the case presented here, the ground-state wave function was obtained by 
the diagonalization of the effective interaction, obtained from the non-local CD Bonn 
NN interaction, while three-body forces were neglected. In the absence of an effective 
theory, which would achieve consistent description of CP-conserving and CP-violating 
observables, the model dependence cannot be completely removed from our results. 
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Figure 10. The polarization contribution to the EDM of ■^He, decomposed into 
isoscalar, isovector and isotensor components. We present contributions for both tt 
and p exchanges; the ^He ground state has been obtained using the CD Bonn NN 
interaction. Convergence with the model-space size for different HO frequencies is 
shown. 



However, the model dependence for the pion exchange is small, since the long-range 
part of the wave function shows little model dependence. Not surprisingly, the short- 
range contributions to the EDM, i.e., p and uj meson exchanges, present a fairly strong 
dependence on the choice of the CP-conserving NN potential model. Nevertheless, 
because ^He is mainly an S'-state nucleus, and the CP-violating Hamiltonian involves 
S- to P-wave transitions, the effect is two fold: (i) the long range contribution, i.e., tt, 
is enhanced, and (ii) the short-range contribution (p and u exchanges) are suppressed. 
Hence, if the CP-coupling constants are of similar magnitude as expected, the vr 
contribution dominates. A quick comparison in Fig. shows that indeed the p (and 
similarly the u) contribution to the ^He EDM is only about 10% of the vr exchange. In 
the hypothesis that the unknown CP-violating coupling constants G'^ are of the same 
order of magnitude for vr, p and u exchanges, we, thus, conclude that the tt contribution 
dominates the EDM, and its value is 

D = (0.024 (5° + 0.023 Gl + 0.027 Gl) e fm. (39) 

This value was obtained after a compilation of all the potential models we have used 
(see Ref. [92]). (Note that in the absence of isospin violation in the Hamiltonian, 
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there is a trivial relationship between the EDM of ^He and that of the triton, i.e., the 
isoscalar and isotensor contributions change sign; the small isospin violation slightly 
breaks this symmetry.) If we compare this result with the neutron dipole moment 
dn = (0.010 G° - 0.010(52) e fm, and deuteron EDM, Ddeut = 0.015 e fm, we see 
immediately that the three EDMs are complementary. Consequently, we can conclude 
that a measurement of these three systems would provide a valuable constraint for the 
theoretical models of CP-violating NN interactions. 

4-2. Charge radii of He isotopes 

Recent advances in the theory of the atomic structure of helium as well as in the 
techniques of isotopic shift measurement made it possible to determine accurately the 
charge radius of ^He [91] and ^He [95]. Precision laser spectroscopy on individual 
^He and ^He atoms confined and cooled in a magneto-optical trap was performed 
and measured the isotope shift between ^He, ^He and ^He. With the help of precise 
quantum-mechanical calculations with relativistic and QED corrections [96] and from 
the knowledge of the charge radius of ^He (1.673(1) [97]), it was possible to determine 
the charge radius of ^He to be 2.054 ± 0.014 fm [9l| and the charge radius of '^He to be 
1.93 ± 3 fm ^95j. The large differences between the ^He and ^He and between the "^He 
and ^He charge radii is due to the extra loosely-bound neutrons in ^He and ^He that 
form a halo [98J . The reduction in charge radius from ^He to ^He indicates a change in 
the correlations of the excess neutrons. 

It is a challenge for ab initio many-body methods to calculate the nuclear radii with 
an accuracy comparable to current experimental accuracy and, thereby, test the nuclear 
Hamiltonians used as the input of ab initio calculations. In Ref. ^0], the ground-state 
properties of ^He, ^He and ^He were calculated within the NCSM using two different 
high-precision NN potentials: the CD-Bonn [2H| and the INOY[29j. 

The ^He calculations were performed both in the Slater determinant basis using 
the Antoine code [57] with model spaces up to A^max = 22 within the two-body effective 
interaction approximation and the Jacobi-coordinate HO basis using the Manyeff 
code [36] with model spaces up to A^max = 20 within either the two-body or the three- 
body effective interaction approximation (with both approximations converging to the 
same result). The ground-state energy and radius convergences are good for both NN 
potentials. The NCSM calculations for ^He and ®He nuclei were performed within 
the two-body effective interaction approximation that allows one to reach much larger 
model-space sizes than within the three-body effective interaction approximation. As 
the radius operator is a long-range operator, it is essential to use as large an HO basis as 
possible. Using the Antoine code, we were able to reach model spaces up to A^max = 16 
and A^rnax = 12 for ^He and ^He, respectively, for a wide range of HO frequencies. 

The point-proton root-mean-square (rms) radii results are summarized in Table [31 
We note that the point-proton rms radius is related to the proton charge rms radius as 
follows [nH (r2) = (r2) - (/?2) _ {Rl)(N/Z), with ((i?2))V2 = 0.895(18) fm [99] (the 
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Table 3. Point-proton (rp) rms radii of ^^^'^He isotopes. The calculated values 
were obtained within the ab initio NCSM [60] . The experimental values are from 
Refs. [nZllMKM]- 
rp [fm] Expt. CD-Bonn 2000 INOY 

^He 1.455(1) 1.45(1) 1.37(1) 

^He 1.912(18) 1.89(4) 1.76(3) 
^He 1.81(3) 1.88(6) 1.74(6) 

charge radius of the proton) and (Rl) = —0.120(5) fm^ [lOOj (the mean-square-charge 
radius of the neutron). In Fig. [TTl we show the model-space size dependence of the rms 
radii for different HO frequencies. A general feature is a decrease of the HO frequency 
dependence with increasing model-space size defined by A^max- In all cases, the rms radii 
exhibit convergence. The ^He point-proton rms radius experimental value is shown as 
a dashed line in top panels of Fig. [HI with the dotted lines indicating the experimental 
error. The CD-Bonn 2000 ^He point-proton rms radius stabilizes at A^max = 16 for the 
HO frequencies of HQ = 9 and 10 MeV, while it is still decreasing for HQ = 8 MeV 
but is increasing for the HO frequencies higher than hfl = 10 MeV. Clearly, the stable 
result is very close to the experimental value. The error was estimated from the HO 
frequency dependence of the A^max = 16 calculations. Based on these results, we arrive 
at the CD-Bonn 2000 point-proton rms radius of 1.89(4) fm that, taking into account 
the error bars, agrees with the experimental value of 1.912(18) fm. We observe a better 
convergence for the INOY NN potential not only for the binding energies, discussed in 
detail in Ref. [60], but also for the radii. This is apparent from Fig. [TTl For this NN 
potential, we find the ^He point-proton rms radius to be 1.76(3) fm. This is significantly 
less than in experiment. Clearly, the INOY NN potential underpredicts both the ^He 
and ^He point-proton rms radii. 

The NCSM ^He point-proton rms radius results are shown in bottom panels of 
Fig. [m for the CD-Bonn 2000 and INOY potentials. From the basis size and the HO 
frequency dependence, a prediction was made in Ref. [60] that the ^He point-proton rms 
radius is 1.88(6) fm based on the NCSM CD-Bonn resuhs. The INOY NN potential 
gives a smaller value, 1.74(6) fm, consistent with the smaller ^He and ^He results. In 
both cases, the ^He point-proton radius is slightly smaller then the corresponding one 
in ^He. The subsequent experimental measurement [95] found the ^He radius of 1.93(3) 
fm, which translates to a point-proton rms radius of 1.81(3) fm, which is slightly lower 
than the NCSM CD-Bonn result but still consistent with it, taking into account the 
theoretical and experimental error bars. 

4-3. Electromagnetic moments of Li and Be isotopes 

Recent developments of both experimental and theoretical techniques have allowed for 
very precise measurements of charge radii and ground-state electromagnetic moments of 
exotic isotopes |101l I102| 11031 1104] . Electric quadrupole and magnetic dipole moments 
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Figure 11. Dependence of the ^He (top panels) and ^He (bottom panels) point-proton 
rms radius on the model-space size for different HO frequencies, obtained using the CD- 
Bonn 2000 (left panels) and the INOY (right panels) NN potentials. The experimental 
values are from Rcfs. [94] and [95] , 



in particular can be determined using an experimental method that is based on the 
nuclear magnetic resonance technique [1031 1104] . These observables reflect, in different 
ways, the evolving nuclear structure along the isotopic chains. 

In response to these and forthcoming experimental programs, extensive ab initio 
NCSM calculations were performed for charge radii and electromagnetic moments of Li 
and Be isostopes [58] in a similar way as for the charge radii of He isotopes discussed 
in the previous subsection. In particular, the CD-Bonn and the INOY NN potentials 
were employed. The calculations were performed up to very large model spaces in a 
wide range of HO frequencies. Efforts were made to quantify the rates of convergence of 
observables. In order to maximize the size of the basis, a two-body effective interaction 
approximation was employed. 

Detailed results for the ground-state energies, charge radii, quadrupole and 
magnetic moments of ^''^''^'^'^^Li and '''^'I'^'i^Be were published in Ref. [58]. Here we 
compare in Fig. [12] the calculated ground-state energies of Li isotopes and magnetic 
moments of Li and Be isotopes obtained using the INOY NN potential with experimental 
trends. As already discussed in Subsection 12. 5[ convergence of eigenenergies with this 
potential is very good and a reliable exponential extrapolation, i.e. using ii^(A^max) = 
Eao + a exp(— feA^inax), can be performed. The ground-state energies of the Li isotopes are 
nicely reproduced by the INOY potential except for ^^Li where lower A^max truncation is 
used compared to lighter Li isotopes as dimensions grow steeply with A. Consequently, 
the ground-state energy extrapolation is more uncertain. Overall, with the exception 
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Figure 12. The Li isotope ground-state energies (left panel) and Li and Be isotope 
magnetic moments (right panel) obtained using the INOY NN potential. The magnetic 
moment experimental values are from Refs. |1031 11041 11051 1106] . 



of the radius of the ^^Li halo ground-state we find a very good agreement between 
NCSM results and recent experiments. The overall trends of all observables are well 
reproduced. Magnetic dipole moments characterized by good convergence properties 
with the NCSM are found in agreement with the experimental trend. Another success 
is the tiny quadrupole moment of ^Li that is known to pose a difficult task for most 
theoretical calculations. In particular, the general failure of three-body models for this 
observable has been blamed on missing antisymmetrization of the valence nucleons and 
the nucleons in the alpha-core. The NCSM correctly reproduces the very small value, 
but with CD-Bonn and INOY giving different signs. Simultaneously, the trend for the 
much larger moments of A = 7 — 11 is nicely reproduced. We note that the ratio 
Q{^^Li)/Q{^Li) is found to be very close to unity, as confirmed recently by very precise 
experimental data |104] . This finding is obtained without a very accurate description 
of the dilute halo structure of ^^Li; a structural feature that we find would require an 
extension of the HO basis used in the standard NCSM. Still, the decrease of the charge 
radius of A = 6 — 9 isotopes is reproduced, although the INOY interaction gives too 
high nuclear densities. 

4.4- Natural and unnatural parity states of^Be and ^^Be 

It is a challenge for nuclear theory to describe odd-A beryllium isotopes from first 
principles. The A = 11 isobar is of particular interest in this respect since it exhibits 
some anomalous features that are not easily explained in a simple shell-model framework. 
Most importantly, the parity-inverted 1/2"*" ground state of ^^Be was noticed by Talmi 
and Unna [107] already in the early 1960s, and it still remains one of the best examples 
of the disappearance of the N = 8 magic number. 

Large-basis ab initio NCSM calculations for ^Be and ^^Be were reported in 
Ref. [108] . Calculations were performed for both natural-parity and unnatural-parity 
states in model spaces up to N^^^ = 9 using four different accurate NN potentials: CD- 
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Figure 13. Excitation spectrum for ^Be (left panel) and ^^Be (right panel) calculated 
using the AV8' interaction in Nj^ax = — 9 model spaces with a fixed HO frequency 
of = 12 MeV. The experimental values are from Refs. [106] and |109| . respectively. 
The AV8' ^Be results obtained by the GFMC method [TUI are shown for comparison. 



Bonn 2000, Argonne V8' [I0l[27], INOY and chiral EFT N^LO. To maximize the model- 
space size, the NNN forces were not included and the two-body effective interaction 
approximation was used. The investigation included a ^^B A^max = 9 calculation with a 
basis dimension of 1.1 x 10^, the largest NCSM diagonalization at that time. 

General features of the excitation energy results are represented in Fig. [131 We 
observe a very reasonable agreement with experimental levels of natural parity, while 
the unnatural-parity states are consistently high in excitation energy. For both parities, 
there is a general trend of convergence with increasing model space. It is clear, however, 
that the relative position of the negative- versus positive-parity states is still not 
converged. Furthermore, when studying the AV8' convergence pattern, it seems as 
if this interaction will predict the positive-parity states at too high excitation energies, 
in particular for ^^Be, even when the calculations will be converged. For ^Be, it was 
found that calculations with the AV8' and N'^LO interactions predict the first-excited 
negative-parity state to be a l/2~, while experiments show that it is a 5/2~. This level 
reversal was also found in the GFMC calculations using AV8'. The INOY interaction, on 
the other hand, gives the correct level ordering, but instead overpredicts the spin-orbit 
splitting. 

The experimental ground state of ^^Be is an intruder 1/2"*" level, while the first 
p-shell state is a l/2~ situated at = 320 keV. The neutron separation energy is 
only 503 keV, and there are no additional bound states. This level-ordering anomaly 
constitutes the famous parity-inversion problem. An important topic of this work has 
been the investigation of the parity inversion found in ^^Be. The large-scale ab initio 
NCSM calculations of Ref. |108j did not reproduce the anomalous 1/2"*" ground state, 
but did observe a dramatic drop of the positive-parity excitation energies with increasing 
model space. Furthermore, the behavior of the INOY results suggested that a realistic 
NNN force might have an important influence on the parity inversion. However, as we 
discuss in the last part of this review, one cannot draw any conclusions without first 
including the extended n—^^Be configurations. 

This conclusion is further re-inforced by the El transition calculations for the strong 
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Table 4. B(E1) values (in [e^fm^]) for the strong ground-state transitions in ^Bc and 
^^Be. The NCSM calculations were performed in the A^max = 8(9) model space for 
negative-(positive-)parity states using the AV8' interaction. Experimental values are 
from [TMfTggifTTUlfTTTj . 



Expt. 0.061(25) 0.0100(84) 0.116(12) 

NCSM 0.033 0.0057 0.0065 



ground-state transitions in ^Be and ^^Be, reported in Ref. [108j and summarized in 
Table |H The strength of the electric dipole transition between the two bound states in 
^^Be is of fundamental importance. This is an observable, which has attracted much 
attention, since it was first measured in 1971 [112] . and again in 1983 [113] . The cited 
value of 0.36 W.u. is still the strongest known transition between low-lying states, 
and it has been attributed to the halo character of the bound-state wave functions. 
Unfortunately, by working in a HO basis, we suffer from an incorrect description of the 
long-range asymptotics, and we would need an extremely large number of basis states 
in order to reproduce the correct form. This shortcoming of the HO basis is illustrated 
by the fact that we obtain a value for the El strength which is 20 times too small. 
When studying the dependence of this value on the size of the model space, we observe 
an almost linear increase, indicating that our result is far from converged. A similar 
increase is observed for ^Be. However, for this nucleus we note that, in the largest model 
space, our calculated El strength is only off by a factor of two compared to experiment. 

p. _|_ Q 

In addition, a consistent result is found for the much weaker ^ — > ^ El transition in 
^Be, where we also obtain a factor of two smaller B{E1) than experiment. These results 
accentuates the anomalous strength observed for ^^Be. As argued in the last part of this 
review, a proper asymptotic behavior of the n—^^Be S-waye configurations that cannot 
be achieved withou extending the NCSM basis by cluster configurations is critical for 
explanation of the B(E1) strength between the two bound states of ^^Be. 

5. Extension to heavier nuclei and larger model spaces 

No-core shell model calculations become computationally intractable for heavier nuclei. 
At least A^max = 6(8) model spaces are required to obtain stable excitation energies 
of the lowest states in the three-body (two-body) cluster approximation (i.e. using 
the three-body (two-body) effective interaction). Still larger spaces are needed for a 
full convergence of binding energies. The m-scheme dimension grows exponentially 
with the number of nucleons A and the truncation level A^max- For this limits 
the presently tractable model space to A'max = 8 [114] . corresponding to an effective 
m-scheme dimension of 6 x 10^. For ^°Ca the dimension of the 8hQ model space is 
2 X 10^^ — well beyond the capabilities of current shell model codes. In order to extend 
applicability of the ab initio NCSM to heavier nuclei and larger model spaces, we must 
resort to approximation schemes. In this section we discuss two such schemes that were 
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developed recently. 



5.1. Importance truncated no-core shell model 

It turns out that many of the m-scheme basis states used in the NCSM calculations 
are irrelevant for the description of any particular eigenstate, e.g., the ground state. 
Therefore, if one were able to identify the important basis states beforehand, one 
could reduce the dimension of the matrix eigenvalue problem without losing predictive 
power. This can be done using an importance truncation scheme based on many-body 
perturbation theory [115] . 

The general concept of the importance truncation is as follows: One starts with a 
reference state I^E'ref), which in the simplest case of closed-shell nuclei could be a single 
harmonic-oscillator Slater-determinant, but, in a general reference state of any 

complexity can be used. Also, excited states can be included, yielding a set of reference 
states. Starting from l^'ref), one can build a many-body space by generating all possible 
n-particle-n-hole (npnh) excitations up to the excitation energy N^g^^f^Q. By increasing 
n (< A), we eventually recover the translationally invariant N^s^^hfl model-space of 
the NCSM. We now estimate the contribution of a given basis state \^^) to the exact 
eigenstate after the diagonalization via many-body perturbation theory. In first-order 
the amplitude of the state is given by 

(40) 

Ci/ ^ref 

where H' is the Hamiltonian of the perturbation and e^, e^ef are the unperturbed 
energies of the two configurations. In the case of closed-shell nuclei with a single- Slater- 
determinant reference state, the unperturbed Hamiltonian is just the one-body harmonic 
oscillator Hamiltonian, Hq = Huo, given, e.g., as the first term on the right-hand side 
of Eq. ([2]), with the Slater determinants 1$^) as eigenstates, i.e., Huol^^) = ^u\^u) and 
-f^Hol^ref) = ^ref l^ref)- The perturbation is given by H' = — ifnO) with the internal 
Hamiltonian Ha given by Eq. ([1]). 

If we restrict ourselves to two-body interactions in Ha ([I]), then H' contains only 
one- and two-body terms such that K,y vanishes for 3p3h and higher-order configurations, 
when starting with a OpOh reference state. In principle, higher orders of perturbation 
theory are required to directly generate states beyond the 2p2h level. Since this becomes 
computationally inefficient, one can resort to an iterative scheme. In a first iteration 
we generate Iplh and 2p2h, states starting from a Slater determinant as the reference 
state; retain those with an importance weight > /tmin! and solve the eigenvalue 
problem in this space. In the next iteration, the dominant components of the ground 
state |\E'o) = l^j, C'l/I^'y) obtained from the diagonalization in the previous step are 
used as reference state, i.e., \^rei) = ^J^'"" C^\^^) with Cj-ef ~ 0.0005. Since this 
state already contains up to 2p2h admixtures, one obtains nonvanishing importance 
weights pUj) for states up to the 4p4h level. After applying the importance truncation 

> Kmin, we solve the eigenvalue problem in the extended space. This cycle can be 




Figure 14. Ground-state energy dependence on the model-space size iVmax for O, 
obtained within the importance-truncated NCSM, using the SRG-N'^LO NN potential 
of Ref. [nil with a cutoff of 2.66 fm~ ^ . The HO frequency hn = 24 MeV was employed. 
The calculation is variational. No NCSM effective interaction was used. The full NCSM 
results were obtained with the code Antoine [57] . 



repeated until the full N^^ax^^ model-space is generated in the limit Kmin = 0. 

Alternative schemes to the above can be developed [116] . For example, one can 
start from a more complex reference state (or a set of reference states) obtained in an 
A^max^^ space. Then by application of Eq. fHUl) with an appropriate re-definition of Hq 
and H', we generate a new basis with components up to 2p2h above the highest npnh 
configurations present in the reference state (s). We can retain only the configurations 
up to, e.g., (A^max + 2)hQ, diagonalize the Hamiltonian in this importance-truncated 
basis, and obtain a new reference state (or set of reference states). This procedure can 
be repeated up to a desired N^^^^hfl space, optimally until convergence is reached. 

To remove a dependence of the results on the cutoff of calculations can 

be performed with Kmin varied. It then becomes feasible to extrapolate to the desired 
i^min = case |115] . 

In Fig. [HI we show an example of ^^O ground-state convergence obtained 
by the application of the importance-truncated NCSM using the above iterative 
procedure |117j . More results for ^He, ^^O and '^°Ca were presented in Ref. |115j . 

It should be noted that the reduction of the basis dimension that can be achieved by 
this method is quite dramatic. For example, importance-truncated NCSM calculations 
for ^^O in the Nj^.^^ = 16 model space reach dimensions of the order of 10^, which can be 
managed easily. This is one order of magnitude less than the dimension of full A^max = 8 
NCSM calculations for the same nucleus. 

The importance-truncated NCSM calculations are completely variational, i.e., they 
provide an upper bound for the ground-state energy of the system, if no NCSM model- 
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space effective interaction is used, as is done in tlie calculations with the soft potentials 
used in Fig. [H] and also in Ref. |115j . If an NCSM effective interaction is used for a 
specific A^max model space, the importance-truncated NCSM calculation is variational 
only with respect to the full NCSM calculation for the same A^max- 

The ab initio NCSM preserves translational invariance, when working with a Slater 
determinant basis, because it uses the HO single-particle states and truncates the model 
space according to the total HO excitation energy Nj^.^^hQ. This is important not only 
for obtaining proper binding or excitation energies, but also for a correct extraction 
of physical wavefunctions, which is crucial, e.g., in the ab initio reaction calculations, 
discussed later in the paper. The importance-truncated NCSM violates the translational 
invariance in a minimal way that can be monitored by looking at the invariance of 
intrinsic properties (e.g., energy, radii, density jTl9j ) with the value of /? in a Lawson 
projection term ( discussed, e.g., after Eq. ( 12^ ). The variational nature and the minimal 
violation of the translational invariance distinguish the importance-truncated NCSM 
from other truncation schemes, such as those used in the coupled-cluster method. This 
was pointed out in the discussion [120[ 1121] , following the publication of Ref. jll5j . 

5.2. Ab initio shell model with a core 

The two-body and three-body cluster approximations presented in subsections 12.3.21 and 
12.3.31 can be generalized further for a larger size of clusters [17 t I122j . For example, if we 
take a six-nucleon system, we may construct a six-body effective Hamiltonian that takes 
into account the full-space six-nucleon correlations. Assuming that it is possible to solve 
the six-body problem in a very large space using two- or three-body interactions, and to 
achieve reasonable eigenenergy convergence, we may treat the obtained solutions as the 
full space results. In fact, the NCSM calculation for the A = 6 system in the iVmax = 12 
space yields nearly converged energies for the lowest states dominated by the = 
components. 

Moreover, if it is possible to solve the six-body problem for A = 6, then it is 
possible to solve the six-body problem for arbitrary A, using the corresponding effective 
Hamiltonian H^^'""''^ obtained in the two-body cluster approximation. This means that 
we can determine for any A-hodj system the secondary effective six-body Hamiltonian 
which accounts for six-body dynamics in the large A^max = 12 space. The six-body 
effective Hamiltonian can be constructed for an arbitrary Pi space with A?"! max < A^max; 
however, the simplest and most instructive case is A^i,max = (single p-shell), i.e., 
when the secondary six-body effective Hamiltonian has the dimension of the two-body 
Hamiltonian in the p-shell. Thus, we can construct an effective two-body Hamiltonian 
for the p-shell, which also contains information about 3-,4-, 5- and 6-body correlations in 
the large A^max = 12 space for the A-hodj system. Because the Pi space has A^'i^max = 0, 
the projection into this space "freezes" four of the A nucleons into fixed single-particle 
configurations, which can be thought of as the " inert core" states in the Standard Shell 
Model (SSM) approach. Consequently, we are dealing with a two-body valence cluster 
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approximation for A — 4 valence nucleons in the p-shell. Remember that we have solved 
the six-body problem for 6 < A < 16. 

To calculate the matrix elements of the secondary effective Hamiltonian, Ti^^, for 
two-body HO basis states, \ap^), we employ the transformation given by Eq. (120|) . where 
the two-body eigenstates \k) and eigenvalues Ek have to be replaced with six-body 
eigenstates and eigenvalues, respectively. 

In the general case of a doubly magic closed shell with two extra nucleons i.e., A = 
6, 18,42, etc., the dimension of the effective Hamiltonian is a 2-body Hamiltonian 
in the p— , sd—, p/— spaces, etc., respectively. This means that the secondary effective 
Hamiltonian does not contain 3- and higher-body terms, even after the exact A-hodj 
cluster transformation. This effective Hamiltonian, which now contains the correlation 
energy of all v4 nucleons, is the correct two-body Hamiltonian to use in a SSM calculation 
with an inert core. The A^ = A — 2 nucleon-spectators fully occupy the shells below 
the valence shell and the total A-hodj wave-function can be exactly factorized as the 
Ac-body "core" and the valence 2-body wave functions. This considerably simplifies the 
calculation of the effective Hamiltonian, because only the OhQ part (Pi-space part) of 
the complete N^g^^hfl wave function needs to be specified. 

Utilizing the approach outlined above, we have calculated effective p-shell 
Hamiltonians for ^Li, using the 6-body Hamiltonians with A^max = 2, 4, .., 12 and f2 = 14 
MeV, constructed from the INOY interaction [22]. The corresponding excitation energies 
of p-shell dominated states and the binding energy of ^Li are shown as a function of 
A^max on the left of Figj2] discussed already in Sect. 12.51 

In the SSM an effective two-body Hamiltonian for a nucleus with mass number A 
is represented in terms of three components: 

ifs^M = ^0 + ^1 + (41) 

where Hq is the inert core part associated with the interaction of the nucleons occupying 
closed shells. Hi is the one-body part corresponding to the interaction of valence 
nucleons with core nucleons, and V^^ is the two-body part referring to the interaction 
between valence particles. It is usually assumed that the core and one-body parts are 
constant for an arbitrary number of valence particles and that only the two-body part 
may contain mass dependence that includes effects of three-body and higher-body 
interactions. 

To represent the Hamiltonian in the SSM format, we develop a valence cluster 
expansion (VCE), 

k=2 

where the lower index, k, stands for the /c-body interaction in the Ov-body valence 
cluster; the first upper index A for the mass dependence; and the second upper index, 
fc + 4, for the number of particles contributing to the corresponding fc-body part. Thus, 
we consider the more general case of allowing the core {k = 0), one-body {k = 1) and 
other /c-body parts to vary with the mass number A. This appears necessary to include 
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the A-dependence of the excitations of the four core nucleons treated in the original 
^max basis space. 

For the A = 6 case the two-body valence cluster (2BVC) approximation is exact: 

where the core part, Hq"^ = is defined as the ground state J'^ = 0"*" energy of ^He 

calculated in the iVmax^^ space with the TBMEs of the primary effective Hamiltonian, 
^iVmax,!^ for A = Q. Then the one-body part, if^'^,is determined as 

Hi'' = n'J - H^. (44) 
The TBMEs of the one-body part, i^f'^, 

{ah\Hl'''\cd)jT = {ea + et)5a,c5b,d (45) 
may be represented in terms of single particle energies (SPE) , ea'- 

el = Ei'U^ja) - Hl'\ el = i?(^He, j„) - H','\ (46) 

where the index a (as well as b, c, and d) denotes the set of single particle HO quantum 
numbers {na,la,ja), the upper index stands for proton (p) and neutron (n), and the 
^(^Li, J) and ^(^He, J) are the NCSM energies of the lowest J[ = 3/2^ and J[ = 1/2^ 
states calculated in the A^max^^ space for the 5-body system using the TBMEs of the 
A = 6 effective Hamiltonian, H^^^f' ^ which includes Coulomb energy. Finally, the 
two-body part ^^'^ is obtained by subtracting the two Hamiltonians: 

v^'' = n'^^-n%. (47) 

It is worth noting that since the Coulomb energy is included in the original Hamiltonian, 
the proton-proton (pp), neutron-neutron (nn) and proton-neutron (pn) T = 1 TBMEs 
of the two-body part, V^^'^, have different values. 

The VCE given by the Eq. fl42l) would require a three-body part V3' of the p-shell 

7 3 

effective interaction Ti^^ to reproduce exactly the NCSM results for A = 7 nuclei: 

Therefore, it is worth knowing how good the 2BVC approximation is for A = 7, as 
well as for A > 7. To test the 2BVC approximation, we have constructed the Ti^'^'^ 
Hamiltonian, using an expansion in terms of zero-, one- and two-body valence clusters, 
i.e., omitting the three-body part: 

nt=''' = H'/ + Hi'' + V^'\ (49) 

In other words, we have first performed NCSM calculations for the k + 4-body systems 
{k = 0, 1, 2) with the H^^f^Y Hamiltonian. Thus, H^'^ is the ^He ' 'core" energy and 
hI'' is the one-body part determined, as in Eqs. fH^ - fj46|) . but with A = 7] and ¥2'^ 
is obtained by subtracting H'^''^ + Hi'' from Tif^'^''^. The resulting parts of the Tif^'^''^ 
Hamiltonian are given in Table [5l 
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Table 5. The core energy, SPEs and pn TBMEs of the effective p-sheh Hamiltonian 
TC^g for the A — 6 and the A — 7 systems in the 2BVC approximation. 







Eq fH3l) 


Ea fl49l) 


Ea fl50l) 






A=6 


A=7 


A=7 




Core energy 


-54.830 


-63.336 


-30.500 




e(P3/2) 


13.922 


10.637 


2.784 




e(Pi/2) 


17.964 


16.355 


5.355 


(P3/2 


^2^''b^/2)^=3,T=0 
^2^*'^ 1^3/2) J=2,T=1 
^2^'^bv2)j=0,T=l 


-2.181 


-2.457 


-19.586 


(P3/2 


2.094 


2.875 


-14.245 


(P3/2 


-2.823 


-3.104 


-3.104 



Comparing several TBMEs for A = 6 and A = 7 (Table [5]), we find that they differ 
considerably. There is a big change separately for the core and one-body parts, but 
weaker changes for the two-body parts, which tend to become larger in magnitude with 
increasing A. We have then performed SSM calculations for the ground state energy of 
^Li, using the zero-, one- and two-body parts in Eq. fl49l) . Namely, the one- and two-body 
parts were employed in a SSM calculation of the ground and excited-state energies of 
the valence nucleons in the p-shell, i.e., OhO, space, to which the ^He core energy, Hq'"^, 
was added, in order to yield the total energies. These calculations were repeated for 
iVmax = 0,2, ...10. Next we carried out NCSM calculations for ^Li with H^ZT^f for 
the same values of A^max- The SSM and NCSM results for the ground-state energy are 
shown in FigJTSl 



> 



0> 



28 
32 
■36 
■40 
■44 
48 



a) ^Li INOY hQ=14 MeV 



Exp 



6 

Nn 



8 10 12 



12 
10 



> 



X A 



b) 


■'Li 


INOY 


hQ=14 


MeV 


2 J 


























• • 




7- »- 














■ -t — 
=1= 




— • 



8 10 12 



Figure 15. a) The ground-state energy, iJgs, of ''Li as a function of iVmax for the 
INOY interaction. The NCSM results with the H^^j'^f^ Hamiltonian are shown by 

' 7 2 

filled circles connected with the solid line. The SSM results with the effective 
Hamiltonian decomposed according to Ea. (|49| are shown by filled circles connected 
with a dashed line. The SSM results with the effective "H^^ Hamiltonian decomposed 
according to Ea. ((50| are shown by filled squares connected with a dashed line, b) 
NCSM (sohd line) and SSM (using Eq.dM]), dashed line) spectra for "^Li. The states 
with spin J are marked by 2 J. 
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It is also of interest to find out what would be the result if we take the fixed core 
and one-body parts at values, which are appropriate for the A = 4 and A = 5 systems, 
respectively, because this is analogous to what is done in the SSM to determine energies 
relative to an inert core. To do this we adopt an alternative two-body VCE, which 
assumes that the core and one-body parts are A- independent, i.e., 

= Ht^' + H'/ + Wt^\ (50) 

similar to the SSM convention given by Eq. fj^TI) . We have then performed another set 
of SSM calculations for A = 7 in the same manner as described previously, but using 
the decomposition given in Eq. (!50!) . To distinguish between the two-body parts of the 
VCE given by the Eqs. (l42l) and (jSOll . we have introduced the new notation, W.^'^, in 
Eq.( l50l) . The Hamiltonian TiJ^ expanded according to the Eq. (|50ll is shown in last 
column of Table [5] and the corresponding results are depicted in FigJTSk) by the solid 
squares connected with a dashed line. Figure [15] indicates that for light systems a 
realistic balance of core, one-body and two-body parts of the effective interaction may 
be achieved only when both the core and one-body parts are calculated according to 
the prescription given above, i.e., Eq.fHOl). Adoption of the decomposition procedure 
with A-independent core and one-body parts leads to a very strong diagonal two- 
body part for the valence nucleons and, subsequently, to drastic overbinding. It is 
obvious, that, in order to compensate for such an effect, one would need to introduce a 
strongly repulsive three-body effective interaction with an unrealistic strength of about 
10 MeV. The VCE with the A-dependent core and one-body parts also yields better 
agreement with the exact NCSM results for the excited states. The corresponding low- 
energy spectra of '^Li obtained with the NCSM and the A-dependent SSM are shown in 
FiglTSb). The differences observed in FiglTSk) and b) for the ground state and excited 
states, respectively, may be attributed to the neglected three-body part of the effective 
interaction at the two-body valence cluster level. 

We have generalized the 2B VC expansion procedure of Eq. ( l49i) for arbitrary mass 
number A, 

and applied it to the A = 7,8,9, and 10 isobars for A^max = 6 in |122j . We have 
found that the three-body and higher-body correlations become more important with 
increasing mass number. There is also a very strong isospin dependence of the obtained 
results. For the highest isospin values the SSM systematically underbinds nuclei in 
comparison to the NCSM and higher-body correlations appear to be small for systems 
containing only valence neutrons. However, there is an opposite effect in the vicinity of 
the N = Z line, where SSM yields considerably more binding energy than the NCSM. 

The analysis of the A = 7 systems allowed us to derive an effective three-body 
Hamiltonian for the p-shell and to give an idea about the strength of the effective three- 
body interaction |122j . 
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6. Cluster overlap functions and spectroscopic factors 

In the ab initio NCSM calculations, we are able to obtain wave functions of low-lying 
states of light nuclei in large model spaces. An interesting and important question is, 
what is the cluster structure of these wave functions. That is we want to understand, how 
much, e.g., an ^Li eigenstate looks like ^He plus deuteron, an ^Be eigenstate looks like 
^He plus ^He, an eigenstate looks like ''Be plus proton, and so on. This information 
is important for the description of low-energy nuclear reactions. To gain insight, 
one introduces channel-cluster form factors (or overlap integrals, overlap functions). 
The formalism for calculating the channel-cluster form factors from the NCSM wave 
functions was developed in Ref. [123j . Here we just briefly repeat a part of the formalism 
relevant to the simplest case, when the lighter of the two clusters is a single nucleon. 

We consider a composite system of A nucleons, e.g., ®B, a nucleon projectile, e.g., a 
proton, and an A — 1-nucleon target, e.g., ^Be. Both nuclei are assumed to be described 
by eigenstates of the NCSM effective Hamiltonians expanded in the HO basis with 
identical HO frequency and the same (for the eigenstates of the same parity) or differing 
by one unit of the HO excitation (for the eigenstates of opposite parity) deflnitions of the 
model space. The target and the composite system is assumed to be described by wave 
functions expanded in a Slater determinant single-particle HO basis (that is obtained 
from a calculation using a shell-model code like Antoine). 

Let us introduce a projectile-target wave function 

(e;...eV2r'r|$|;^-;:^)^*'^ = ^(jWiMi|JM)(/m,im,|jm)^ 

X Yi^,{f)XmA^i ■ --U-M - lahM,) , (52) 

where (^i . . . ■^^-2!^ — la/iMi) and Xms are the target and the nucleon wave functions, 
respectively. Here, / is the channel relative orbital angular momentum, C,i are the target 
Jacobi coordinates deflned in Eq. (jlj) and r = (ri + r2 + . . . + r^-i) — ta] describes 
the relative distance between the nucleon and the center of mass of the target. The spin 
and isospin coordinates were omitted for simplicity. 
The channel cluster form factor is then deflned by 

^,.A-..i.(r) = {AXJ\M>\^;^^;6;), (53) 

with A the antisymmetrizer and \AXJ) an eigenstate of the A-nucleon composite 
system (here ^B). It can be calculated from the NCSM eigenstates obtained in the 
Slater-determinant basis from a reduced matrix element of the creation operator. The 
derivation is as follows. First, we use the relation (fT3l) for both the composite A-nucleon 
and the target A — 1-nucleon eigenstate. With the help of HO transformations: 

^{LMlm\Qq)ipNLM{RcM)'^nlm{rA) = 

Mm 

J2 {n'l'N'L'Q I NLnlQ) _j_ {I'm'L'M' \ Qq) 

n'Vm'N'UM' 

X V'n'i'm'(^A-l)v'Af'L'Af'(^o) , (54) 
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we obtain 

SD(AAJ|^<l>|;|-^;^^^;n/)sD = (n/00/|00n//)^ (AAJ|^$j;^-;;^)^; n/) , (55) 

with a general HO bracket due to the CM motion. The nl in (l55ll refers to a replacement 
of Sr by the HO Rni{r) radial wave function. Second, we relate the SD overlap to a 
linear combination of matrix elements of a creation operator between the target and 
the composite eigenstates sD{A^J\(4iij\^ ~ lttA)sD- The subscript SD refers to the fact 
that these states were obtained in the Slater determinant basis. Such matrix elements 
are easily calculated by shell model codes. The result is 

{i^,j);ali ^ (n/00/|00n//)_i_ J 

X sD(^AJ||al,^.p-la/i)sD • (56) 

The eigenstates expanded in the Slater determinant basis contain CM components. A 
general HO bracket, whose value is simply given by 

then appears in Eq. flS^ in order to remove these components. The Rni{r) in Eq. fl5B|) 
is the radial HO wave function with the oscillator length parameter b = ^ ^-i^^o ; where 
m is the nucleon mass. 

A conventional spectroscopic factor is obtained by integrating the square of the 
cluster form factor: 

A generalization for projectiles (= the lighter of the two clusters) with 2, 3 or 
4 nucleons is straightforward, although the expressions become more involved. In all 
cases, the projectile is described by a wave function expanded in the Jacobi coordinate 
HO basis, while the composite and the target eigenstates are expanded in the Slater 
determinant HO basis. Full details are given in Ref. |123j . 

7. S-factors of capture reactions 

The overlap functions introduced in the previous section can be used as a strating point 
for description of low-energy 7-capture reactions important for nuclear astrophysics. 

7.1. '^Be(p,-ffB 

The ^Be(p,7)^B capture reaction serves as an important input for understanding the 
solar neutrino flux |124j . Recent experiments have determined the neutrino flux emitted 
from with a precision of 9% [125] . On the other hand, theoretical predictions have 
uncertainties of the order of 20% |126[ I127j . The theoretical neutrino flux depends on 
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the ^Be(p,7)^B S-factor. Many experimental and theoretical investigations studied this 
reaction. 

In this subsection, we discuss a calculation of the ''Be(p,7)^B S-factor starting 
from ab initio wave functions of ^B and ^Be. It should be noted that the aim of ab 
initio approaches is to predict correctly absolute cross sections (S-factors), not only 
relative cross sections. The full details of our ^Be(p,7)®B investigation were published 
in Refs. [T28lfT29] . 

Our calculations for both ^Be and ^B nuclei were performed using the high-precision 
CD-Bonn 2000 NN potential [28] in model spaces up to lOhQ (A'^max = 10) for a wide 
range of HO frequencies. From the obtained ^B and ''Be wave functions, we calculate 
the channel cluster form factors (overlap functions, overlap integrals) 
as discussed in the previous section. Here, A = 8, / is the channel relative orbital 
angular momentum and r = [-^^ {fi-\-r2 + --- + ^A-i) — ^4] describes the relative 
distance between the proton and the center of mass of ^Be. The two most important 
channels are the waves, / = 1, with the proton in the j = 3/2 and j = 1/2 states, 
j = / + s, s = 1/2. In these channels, we obtain the spectroscopic factors of 0.96 and 
0.10, respectively. The dominant j = 3/2 overlap integral is presented in Fig. [16] by 
the full line. The lOhQ model space and the HO frequency of hfl = 12 MeV were used. 
Despite the fact that a very large basis was employed in the present calculation, it is 
apparent that the overlap function is nearly zero at about 10 fm. This is a consequence 
of the HO basis asymptotic behavior. As already discussed, in the ab initio NCSM, the 
short-range correlations are taken into account by means of the effective interaction. 
The medium-range correlations are then included by using a large, multi-hfl HO basis. 
The long-range behavior is not treated correctly, however. The proton capture on ^Be 
to the weakly bound ground state of ®B associated dominantly by the El radiation is 
a peripheral process. In order to calculate the S-factor of this process we need to go 
beyond the ab initio NCSM, as done up to this point. We expect, however, that the 
interior part of the overlap function is realistic. It is then straightforward to find a 
simple correction to the asymptotic behavior of the overlap functions, which should be 
proportional to the Whittaker function. 

One possibility we explored utihzes solutions of a Woods-Saxon (WS) potential. In 
particular, we performed a least-square fit of a WS potential solution to the interior 
of the NCSM overlap in the range of — 4 fm. The WS potential parameters were 
varied in the fit under the constraint that the experimental separation energy of ^Be-|-p, 
Eq = 0.137 MeV, was reproduced. In this way we obtain a perfect fit to the interior of 
the overlap integral and a correct asymptotic behavior at the same time. The result is 
shown in Fig. [16] by the dashed line. 

Another possibility is a direct matching of logarithmic derivatives of the NCSM 
overlap integral and the Whittaker function: ■^ln{rQij{r)) = ■^ln{CijW-ri,i+i/2{'^kor)), 
where r] is the Sommerfeld parameter, ko = a/2/IEo/^ with /i the reduced mass and Eq 
the separation energy. Since asymptotic normalization constant (ANC) Cij cancels out, 
there is a unique solution at r = Rm- For the discussed overlap presented in Fig. [THl we 
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Figure 16. Overlap function, 
r£|(r), for tlie ground state of 
witli tiie ground state of ^Be plus 
proton. See the text for details. 



Figure 17. The '^Be{p,j)^B S- 
factor from the corrected NCSM 
overlap functions. Experimental 
values are from Refs. [1321 1133j . 



found Rjn = 4.05 fm. The corrected overlap using the Whittaker function matching is 
shown in Fig. [12] by a dotted hne. In general, we observe that the approach using the 
WS fit leads to deviations from the original NCSM overlap starting at a smaller radius. 
In addition, the WS solution fit introduces an intermediate range from about 4 fm to 
about 6 fm, where the corrected overlap deviates from both the original NCSM overlap 
and the Whittaker function. Perhaps, this is a more realistic approach compared to the 
direct Whittaker function matching. In any case, by considering the two alternative 
procedures we are in a better position to estimate uncertainties in our S-factor results. 

In the end, we re-scale the corrected overlap functions to preserve the original NCSM 
spectroscopic factors (Table 2 of Ref. |128j ). In general, we observe a faster convergence 
of the spectroscopic factors than that of the overlap functions. The corrected overlap 
function should represent the infinite space result. By re-scaling a corrected overlap 
function obtained at a finite A^max, we approach faster the infinite space result. At the 
same time, by re-scaling we preserve the spectroscopic factor sum rules. 

The S-factor for the reaction ^Be(p,7)^B also depends on the continuum wave 
function, R^f ■ As the capture reaction calculations were first performed before the 
extension of the NCSM to describe continuum wave functions introduced in Sect. [H the 

(c) 

continuum wave functions Rlj were obtained for s and d waves from a WS potential 
model. Since the largest part of the integrand stays outside the nuclear interior, one 
expects that the continuum wave functions are well-described in this way. In order to 
have the same scattering wave function in all the calculations, we chose a WS potential 
from Ref. |130j that was fitted to reproduce the p-wave 1+ resonance in ^B. It was argued 
|131] that such a potential is also suitable for the description of s- and ci-waves. We 
note that the S-factor is very weakly dependent on the choice of the scattering-state 
potential (using our fitted potential for the scattering state instead changes the S-factor 
by less than 1.5 eV b at 1.6 MeV with no change at MeV). 
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Figure 18. The overlap function, rg(r), for tlie ground state of ^Be witli '^He plus a. 
See the text for details. 

Our obtained S-factor is presented in Fig. [TTl where the contribution from the two 
partial waves are shown together with the total result. It is interesting to note a good 
agreement of our calculated S-factor with the recent Seattle direct measurement |132j . 

In order to judge the convergence of our S-factor calculation, we performed a 
detailed investigation of the model-space-size and the HO frequency dependencies. We 
used the HO frequencies in the range from HQ = 11 MeV to hQ = 15 MeV and the 
model spaces from 6hQ to lOhQ. By analysing these results, we arrived at the S-factor 
value of ^17(10 keV) = 22.1 ± 1.0 eV b. 

7.2. ^He(a,-ifBe and ^ H (a f Li 

The ^He(a,7)'^Be capture reaction cross section has been identified as the most 
important uncertainty in the solar model predictions of the neutrino fluxes in the p- 
p chain [127]. We investigated the bound states of ''Be, ^He and "^He within the ab initio 
NCSM and calculated the overlap functions of ^Be bound states with the ground states 
of ^He plus ^He as a function of separation between the ^He and the a particle. The 
obtained p-wave overlap function of the ^Be 3/2~ ground state is presented in Fig. [18] 
by the full line. The dashed lines show the corrected overlap function obtained by the 
least-square fits of the WS parameters done in the same way as in the ^B<-i>''Be-|-p case. 
The corresponding NCSM spectroscopic factors obtained using the CD-Bonn 2000 in 
the lOhfl model space for ^Be {12hQ for ^'"^He) and HO frequency of HQ = 13 MeV are 
0.93 and 0.91 for the ground state and the first excited state of ^Be, respectively. We 
note that contrary to the ^B^'^Be-|-p case, the "^Be^^He-l-a; p-wave overlap functions 
have a node. 

Using the corrected overlap functions and a ^He-|-a scattering state obtained using 
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Figure 19. The full line shows the ^He(a,7)'^Be (left) and ^H(a,7)'^Li (right) S- 
factors obtained using the NCSM overlap functions with corrected asymptotics. The 
dashed lines show the ^Be (left) and ''Li (right) ground- and the first excited state 
contributions. The calculations were done using the CD-Bonn 2000 NN potential and 
the Whn model space for ^Be and "^Li {Uhn for ^H, ^He and ^He) with the HO 
frequency of hfl ~ 13 MeV. 



the potential model of Ref. |134j . we calculated the ^He(a,7)^Be S-factor. Our lOhfl 
result is presented in the left panel of Fig. [191 We show the total S-factor as well as the 
contributions from the capture to the ground state and the first excited state of ^Be. 
By investigating the model space dependence for 8hQ and lOhQ spaces, we estimate 
the ^He(Q;,7)'^Be S-factor at zero energy to be higher than 0.44 keV b, the value that 
we obtained in the discussed case is shown in the left panel of Fig. [191 Our results are 
similar to those obtained by K. NoUett |135j . using the variational Monte Carlo wave 
functions for the bound states and potential-model wave functions for the scattering 
state. We note that the recent evaluation |136j that took into account precise LUNA 
measurements at low energies [1371 1138] found the S-factor of 0.580 ± 0.054 keV b at 
zero energy. 

An important check on the consistency of the ^He(a,7)^Be S-factor calculation is 
the investigation of the mirror reaction ^H(a,7)''Li, for which more accurate data exist 
|139j . Our results obtained using the CD-Bonn 2000 NN potential are shown on the 
right of Fig. [191 It is apparent that our '^H(a,7)'^Li results are consistent with our 
^He(a,7)^Be calculation. We are on the lower side of the data and find an increase of 
the S-factor as we increase the size of our basis. 

More details on the ab initio NCSM investigation of the ^He(a,7)'^Be and ^H(a,7)^Li 
S-factors are given in Ref. }140] . 

8. Light nuclei as open systems 

Nuclei are open quantum systems with bound states, unbound resonances, and 
scattering states. A realistic ab initio description of light nuclei with predictive 
power must have the capability to describe all the above classes of states within a 
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unified framework. Over the past decade, significant progress has been made in our 
understanding of the properties of the bound states of hght nuclei starting from realistic 
NN and NNN interactions. This progress was also in part due to the development of the 
ab initio NCSM, as described in the previous sections of this review. The solution of the 
nuclear many-body problem is even more complex, when scattering or nuclear reactions 
are considered. Ab initio calculations for scattering processes involving more than four 
nucleons overall are challenging and still a rare exception [19]. Even calculations of 
resonant states are quite complicated |141j . The development of an ab initio theory of 
low-energy nuclear reactions on light nuclei is key to further refining our understanding 
of the fundamental nuclear interactions among the constituent nucleons and providing, 
at the same time, accurate predictions of crucial reaction rates for nuclear astrophysics. 

The use of the harmonic oscillator (HO) basis in the NCSM results in an incorrect 
description of the wave-function asymptotic and a lack of coupling to the continuum. 
The first applications of the NCSM to the calculation of nuclear reactions that were 
described in Sections [7] required a phenomenological correction of the asymptotic 
behavior of the overlap functions. A fully ab initio approach to nuclear reactions based 
on the NCSM requires an extension or a modification of the NCSM basis. 



8.1. Coupling of ab initio no-core shell model with the resonating group method 

The resonating-group method (RGM) is a microscopic cluster technique in which the 
many-body Hilbert space is spanned by cluster wave functions describing a system of two 
or more clusters in relative motion. Here, we will limit our discussion to the two-cluster 
RGM, which is based on binary-cluster channel states of total angular momentum J, 
parity vr, and isospin T, 

|$,^7) = [( \A-aa^I,^^T^) laa^I^T^) f^^ 

X Y^{rA-a,a) — —■ 59 

J rrA-a,a 

In the above expression, ( \ A—aaiI^^Ti) and \aa2l2'^T2) are the internal (antisymmet- 
ric) wave functions of the first and second clusters, containing A— a and a nucleons 
(a<74), respectively. They are characterized by angular momentum quantum numbers 
Ji and I2 coupled together to form channel spin s. For their parity, isospin and addi- 
tional quantum numbers we use, respectively, the notations tTj, Tj, and ««, with i = 1,2. 
The cluster centers of mass are separated by the relative coordinate 

A-a , A 



'^A—a,a '"A— a,a^yl— a,a . / ■ ' 

A — a ^-^ a 

i=l j=A—a+l 

where {r^, i = 1, 2, ■ ■ ■ , A} are the A single-particle coordinates. The channel states ( l59ll 
have relative angular momentum I. It is convenient to group all relevant quantum 
numbers into a cumulative index v = {A—aail^^Ti, aa2-^2^^^2; si}. 
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The former basis states can be used to expand the many-body wave function 
according to 

l^j^T^ = ^ I drr'^^^A. . (61) 



As the basis states fl59|) are not anti-symmetric under exchange of nucleons belonging 
to different clusters, in order to preserve the Pauli principle one has to introduce the 
appropriate inter-cluster anti-symmetrizer, schematically 



p 

Here the sum runs over all possible permutations P that can be carried out 
among nucleons pertaining to different clusters, and p is the number of interchanges 
characterizing them. The coefficients of the expansion are the relative-motion wave 
functions g^'"^{r), which represent the only unknowns of the problem. To determine 
them one has to solve the non-local integro-differential coupled-channel equations 

E / drr^ni:j^{r\r)-EAfnr',r)]'-^ = 0, (63) 

where the two integration kernels, the Hamiltonian kernel, 

ni:j{r\ r) = X,HA. \^t7) , (64) 

and the norm kernel, 

K''7ir',r) = {<l>i:j\ AMu KJ) , (65) 

contain all the nuclear structure and anti-symmetrization properties of the problem. In 
particular, the non-locality of the kernels is a direct consequence of the exchanges of 
nucleons between the clusters. We have used the notation E and H to denote the total 
energy in the center-of-mass frame, and the intrinsic A-nucleon microscopic Hamiltonian, 
respectively. 

The formalism presented above can be combined with the ab initio NCSM as follows. 
First, we note that the Hamiltonian can be written as 

H = r,el(r) + Hel + Voir) + H^^A^a) + H^a) , (66) 

where Hi^a-o) and -f^(a); the {A— a)- and a-nucleon intrinsic Hamiltonians, respectively, 
Tj.ei(r) is the relative kinetic energy and Vrei is the sum of all interactions between 
nucleons belonging to different clusters after subtraction of the average Coulomb 
interaction between them, explicitly singled out in the term Vcir) = Zi^Z2ue'^/r {Zi^ 
and Z2v being the charge numbers of the clusters in channel v). We use identical 
realistic potentials in both the cluster's Hamiltonians and inter-cluster interaction Vrei- 
Accordingly, \A—aaiI^^Ti) and \aa2ipT2) are obtained by diagonalizing H(^A-a) and 
respectively, in the model space spanned by the NCSM basis. If the adopted 
potential generates strong short-range correlations, we derive consistent NCSM effective 
interactions. While the cluster eigenstates are obtained by employing the usual NCSM 
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effective interaction, in place of the bare potential entering Vrei we adopt a modified 
effective interaction, which avoids renormalizations related to the kinetic energy. At the 
two-body cluster level this is given by V^'g^ = H2ea — -f^Ls; where H2es is the effective 
Hamiltonian derived from if^' = Hq2+V(2, with = — mfi^r ^/A, compare to Eq. (123|) . 
Note that Vggg Vn in the limit A^max — > 00. 

As indicated by the presence of the norm kernel A/'^^"^(r', r), Eq. fl63l) does 
not represent a system of multichannel Schrodinger equations, and gf,'"^{r) do not 
represent Schrodinger wave functions. This short-range non-orthogonality, induced by 
the non-identical permutations in the inter-cluster anti-symmetrizers, can be removed 
by introducing normalized Schrodinger wave functions 

J. Z_^y y 

where is the square root of the norm kernel, and applying the inverse-square root 
of the norm kernel, Af'^ ., to both left and right-hand side of the square brackets in 
Eq. fl63|) . This procedure, explaned in more detail in Ref. |143j . leads to the system of 
multichannel Schrodinger equations: 

[T.i(r) + Vc{r) - (E - eS'^^ - E^^"^^ 



dr'r'^ W;^7{r, r') ^^^L^ = q, (68) 



where Ea- ^' is the energy eigenvalue of the i-th cluster {i = 1,2), and W'l^Jir' , r) is the 
overall non-local potential between the two clusters, which depends upon the channel 
of relative motion, while it does not depend upon the energy. 

So far we have fully developed and tested this formalism in the single-nucleon 
projectile basis, i.e., for binary-cluster channel states (1591) with a = 1 (with channel 
index z/ = {A—1 ail^^Ti; 1 si}). In this model space, the norm kernel is rather 
simple and is given by 



A-l 



K''7ir', r) = {<l>i:j\ (69) 

i=l 

= 5.,. ^J^:l_Zll _ _ 1) Vi?„,,,(r')i?„,(0 

r^' ly ^^^^^ 

n'n 

x<$^:j|p^-m|C^>, (70) 

where it is easy to recognize a direct term, in which initial and final state are identical 
(corresponding to diagram (a) of Fig. [20]) . and a many-body correction due to the 
exchange part of the inter-cluster anti-symmetrizer (corresponding to diagram (6) of 
Fig. [20]) . We note that in calculating the matrix elements of the exchange operator 
Pa-\,a we replaced the delta function of eq. fISOl) with its representation in the HO 
model space. This is appropriate whenever (as it is the case of Pa-\,a) the operator is 
short-to-medium range. The presence of the inter-cluster anti-symmetrizer affects also 
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(a) (6) 

Figure 20. Diagrammatic representation of: (a) "direct" and (6) "exchange" 
components of the norm kernel; (c and d) "direct" and (e) "exchange" components of 
the potential kernel. The first group of circled black lines represents the first cluster, 
the bound state of A—1 nucleons. The separate red line represents the second cluster, 
in the specific case a single nucleon. Bottom and upper part of the diagram represent 
initial and final states, respectively. 



the Hamiltonian kernel, and, in particular, the matrix elements of the interaction. For 
a NN potential one obtains a direct term involving interaction and exchange of two 
nucleons only (see diagrams (c) and (d) of Fig. [201) . and an exchange term involving 
three-nucleons. Diagram (e) of Fig. [20] describes this latter term, in which the last 
nucleon is exchanged with one of the nucleons of the first clusters, and interacts with 
yet another nucleon. For more details on the integration kernels in the single-nucleon 
projectile basis we refer the readers to Ref. |143] . 

Being translationally-invariant quantities, the norm and Hamiltonian kernels can 
be "naturally" derived working within the NCSM Jacobi-coordinate basis. However, by 
introducing Slater-determinant channel states of the type 

\A-aaiIiTi)^^ \aa2hT2) )^'''^^ 



\^vn /SD 



X F,(i?, 



{a) \ 
cm. ) 



(71) 



in which the eigenstates of the (A— a)-nucleon fragment are obtained in the SD basis 
(while the second cluster is still a NCSM Jacobi-coordinate eigenstate), it can be easily 
demonstrated that translationally invariant matrix elements can be extracted from those 
calculated in the SD basis of Eq. flTTj) by inverting the following expression: 



E 



jTTrrp 



X 



^£7j2(_1)(^+^-'-^') 



NL 



L J 




X {nJrNLilOOnii)^ {n'/^N Li\00n' i' i') ^ . 



(72) 



Here Ot.i. represents any scalar and parity-conserving translational-invariant operator 
{Ot.i. = A, AHA, etc.). We exploited both Jacobi-coordinate and SD channel states to 
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Table 6. Calculated ''H and ^He g.s. energies (in MeV), n-^H and p-'^He phase shifts 
(in degrees), and n-'^H total cross section (in barns) for increasing iVmax at hil ~ 18 
MeV, obtained using the Mowfe NN potential [SI Ell- The n-^H (p-^He) scattering 
results were obtained in a coupled-channel calculation including only the g.s. of the 
(•^He) nucleus (i.e., the channels u = {3 g.s. l^^^; si}), see text for explanation 

of notation. 
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verify our results. The use of the SD basis is computationally advantageous and allows 
us to explore reactions involving p-shell nuclei. 

8.2. Ah initio many-body calculations of nucleon-nucleus scattering 

The two-cluster NCSM/RGM formalism outlined in the previous section, can be used to 
calculate nucleon-nucleus phase shifts below three-body break threshold, by solving the 
system of multi-channel Schrodinger equations (l68ll with scattering boundary conditions. 
Results for neutrons scattering on ^H, ^He and ^°Be and protons scattering on ^'^He, 
using realistic NN potentials, were presented in Refs. [142] and [143] . In the following 
we review part of these calculations. 

To study the behavior of our approach with respect to the HO model space, we 
performed NCSM/RGM scattering calculations for the A = 4 and 5 systems, using the 
Viowk NN potential [31], which is "soft" and can be treated as "bare", and the N^LO 
NN interaction [30] , which generates strong short-range correlations, thus requiring the 
use of effective interactions. In particular, for these convergence tests, we restricted 
our binary-cluster basis to target-nucleon channel states with the target in its g.s. 
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Figure 21. Calculated n-a phase shifts using the Viowk [511 [Hi] and N^LO [3D] (right 
panel) NN potentials at Hil = 18 and 19 MeV, respectively: (left panel)dependence on 
-^max of the n-a(g.s.) results; (right panel) influence of the lowest six excited states 
(0+0, 0^0, 1-0, 1-1, 2-0, 2-1) of the a particle. 



(corresponding to channel indexes of the type u = {A — 1 g.s. /{"^Ti; 1^ i; si}). A 
sample of the results obtained for Viowfc is presented in Table [6] and in the left panel 
of Fig. [2TJ The overall convergence is quite satisfactory, with a weak dependence on 
^max- The slightly larger differences presented by the phase shifts of small magnitude 
in Table [6] is in part a reflection of the sharp cutoff function used to derive the adopted 
version of V^iowfc from AV18 with a cutoff A = 2.1 fm~^. As emphasized by the left 
panel of Fig. [211 presenting the n-a scattering phase shifts, the convergence rate for 
N^LO (achieved by using two-body effective interactions tailored to the HO model-space 
truncation), is much slower than that obtain with Viowk- However, a gradual suppression 
of the difference between adjacent A/'max values with increasing model-space size is visible, 
although the pattern is somewhat irregular for the P phase shifts. Although not shown, 
the p-a phase shifts present analogous convergence properties. The situation for the 
n-^B. and p-^He systems is similar, and is shown in Fig. [22j Results in the largest model 
spaces (A^max = 17 and 19) are very close. An additional sign of convergence is provided 
by the rather good frequency independence presented by the A^max = 19 n-^B. phase 
shifts. 

The two panels of Fig. [22] show also the results (-1- symbols) obtained for the same 
N^LO NN potential [30] by Deltuva and Fonseca [HI I145j from the solution of the 
Alt, Grassberger and Sandhas (AGS) equations. The discrepancy, increasing with the 
energy, between the two calculations highlights the influence played by closed channels 
not included in our basis states, that is, target-nucleon channel states with the target 
above the /f^ = |^ g.s., and 2+2 configurations, both of which are taken into account in 
the AGS results. Because these states correspond to the breakup of the A = 3 system, 
it is not feasible to include them in the current version of the NCSM/RGM approach, 
which so far has been derived only in the single-nucleon projectile basis. However, we are 
planning on extending our approach to be able to account for the target breakup, and 
these development will be discussed in future publications. For the time being we can 
explore the effects of the virtual excitations of the target on the nucleon-a scattering, 
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Figure 22. Calculated A = A phase shifts using the N^LO NN potential [30] 
compared to AGS results of Refs. p!8l 1145] . In particular: (left panel) n-^H results 
for TVinax = 19 and hn = 19, and 22 MeV; (right panel) p-^Re resuhs for ^ifJ = 22 
IMeV and iVmax = 17, and 19. AU NCSM/RGM results were obtained in a coupled- 
channel calculation including only the g.s. of the ('^He) nucleus (i.e., the channels 
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Figure 23. Calculated phase shifts for (left panel) n-a and (right panel) p-a 
scattering, using the N'^LO NN potential ^U\, compared to an i?-matrix analysis of 
data (+) [146]. Theoretical results include the "He g.s., 0+0, O'O, I'O, I'l, 2-0, and 
2^1 states. 



where much higher energies are needed in order to break up the a particle. This is 
done in the right panel of Fig. [211 which presents the effect of the inclusion of the first 
six excited states of ^He on the and ^^3/2 scattering phase shifts. The use 

of the five different combinations of ground and excited states shown in the legend of 
Fig. [21] indicates that the is well described already by coupled-channel calculations 
with g.s. and first 0"'"0 (the ^5*1/2 phase shifts obtained in the four larger Hilbert spaces 
are omitted for clarity of the figure). On the other hand, the negative-parity excited 
states have relatively large effects on the P phase shifts, and, in particular, the 0~0, 1~0 
and 1^1 states mostly on the ^Pi/2 phase shifts, whereas the 2-0 and 2-1 states on the 
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Figure 24. Calculated n-^'~'Be phase shifts as a function of the relative kinetic energy 
in the cm. frame i?kin, using the CD-Bonn NN potential [28] at = 13 MeV: (left 
panel) and (right panel) ^1^5/2 resuhs. The NCSM/RGM results were obtained 

using n+^°Be configurations with A^max = 6 g.s., 2^, 2^, and 1^ states of ^"Be. The 
obtained scattering length is +10.7 fm. 



The comparison with an accurate i?-matrix analysis of the nucleon-a scatter- 
ing |146j . presented in Fig. [231 reveals that for both neutron (left panel) and proton 
(right panel) projectiles we can describe quite well the ^5*1/2 and, qualitatively, also the 
^D3/2 phase shifts, using the N^LO NN potential. On the other hand, the same interac- 
tion is not able to reproduce well the two P phase shifts, which are both too small and 
too close to each other. This lack of spin-orbit splitting between the and ^Pi/2 

results can be explained by the omission in our treatment of the NNN terms of the chiral 
interaction, which would provide an additional spin-orbit force. This sensitivity of the 
P phases to the strength of the spin-orbit force corroborated by the differences among 
the Viowfc and N^LO results in the left panel of Fig. [211 ^Pi/2 and ^^3/2 are both larger 
and more separated for Viowfc- The good agreement of the N^LO phase shifts with 
their Vowfc analogues and with the i?-matrix analysis can be credited to the repulsive 
action (in this channel) of the Pauli exclusion principle for short nucleon-a distances, 
which has the effect of masking the short-range details of the nuclear interaction. 

Figure [21] highlights one of the promising aspects of the NCSM/RGM approach, 
that is the ability (through the use of SD channel states) to perform ab initio scattering 
calculations for p-shell nuclei. The ^5*1/2 (left panel) and ^-D5/2 (right panel) n-^^Be phase 
shifts were obtained in a A^max = 6, hQ = 13 MeV HO model space. The inclusion of the 
2^ excited state of ^°Be has a significant effect on the S and more importantly on the 
D phase, where it is essential for the appearance of a resonance below 3 MeV. We note 
that a resonance has been observed at ~ 1.8 MeV with a tentative spin assignment of 
(5/2,3/2)+ [109J. The further addition of the 2^ and, especially, the 1^ excited states 
produces rather weak differences. We have also extracted the scattering length for the 
^5'i/2 partial wave and found a result of +10.7 fm, which is comparable to the value 
of +13.6 fm obtained by Descouvemont in Ref. |147] . when fitting the experimental 
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Table 7. Calculated energies (in MeV) of the ^"Be g.s. and of the lowest negative- 
and positive-parity states in ^^Be, obtained using the CD-Bonn NN potential [2S] at 
hfl = 13 MeV. The NCSM/RGM results were obtained using n+^'^Be configurations 
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binding energy of ^^Be. 

8.3. Parity-inverted ground state of^^Be 

Although we mainly described its scattering applications, the NCSM/RGM is a powerful 
tool also for structure calculations, particularly for loosely-bound systems. By imposing 
bound-state boundary conditions to the set of coupled channel Schrodinger equations 
of Eq. ( l68l) . we tested the performance of our single-nucleon projectile NCSM/RGM 
formalism for the description of one-nucleon halo systems. In particular, because of the 
well-known parity-inversion between its two bound states with respect to the predictions 
of the simple shell model |107j . the ^^Be nucleus represents an excellent test ground for 
our approach. 

As discussed in Sec. 14. 4[ large-scale ab initio NCSM calculations with several 
accurate NN potentials of the ^^Be low-lying spectrum were not able to explain its 
g.s. parity inversion |108j . The explanation for these results can be related to two 
main causes: (i) the size of the HO basis was not large enough to reproduce the 
correct asymptotics of the n-^°Be component of the 11-body wave function; and (ii) 
the NNN force, not included in the calculation, plays an important role in the inversion 
mechanism. The second hypothesis was corroborated by the results obtained with the 
INOY NN potential [29J, which produced the lowest excitation energy of the 1/2+ state 
compared to other NN potentials. By studying the ^^Be bound states in a NCSM/RGM 
model space spanned by the n-^°Be channel states with inclusion of the Nmax = 6 g.s., 
2f, 2^, and 1^ of ^°Be, we are now in the position to address the first hypothesis. Indeed, 
the correct asymptotic behavior of the n-^°Be wave functions is described naturally in 
the NCSM/RGM approach. 

The energies of the lowest 1/2+ and 1/2^ states of ^^Be obtained in the NCSM 
and in the NCSM/RGM calculations, using the same CD-Bonn NN interaction [28] at 
HQ = 13 MeV adopted in Ref . [108] , are presented in Table [3 The relatively small 
differences between the A^max = 6/7 and A^max = 8/9 NCSM results, seems to indicate a 
reasonable degree of convergence for these calculations. The l/2~ state appears to be the 
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Table 8. Mean values of the relative kinetic and potential energy and of the internal 
Be energy in the "Be 1/2+ ground state. All energies in MeV. NCSM/RGM 
calculation as in Table [71 See the text for further details. 
NCSM/RGM (r,ei) {W) Ep°Be(g.s.,ex.)] E,,, 

Model Space 16.65 -15.02 -56.66 -55.03 

Full 6.56 -7.39 -57.02 -57.85 



g.s., and the 1/2+ state is about 2.8 MeV above the n-^°Be threshold. A comparison 
to the NCSM/RGM calculations (obtained in a model space including g.s., 2f, 2^, 
and li states of ^°Be) shows a rough agreement for the 1/2" state, whereas for the 
1/2+ state one observes a dramatic difference (~3.5 MeV) in the energy. The l/2~ 
and 1/2+ NCSM/RGM states are both bound and the 1/2+ state is the g.s. of "Be. 
Correspondingly, we obtain a B(E1; ^ — ^ I"*") value of 0.18 fm^, which is not far 
from experiment. 

To understand the mechanism that makes the 1/2+ state bound in the 
NCSM/RGM, we evaluated mean values of the relative kinetic and potential energies 
as well as the mean value of the ^''Be energy, and compared them to those obtained by 
restricting all the integration kernels within the HO model space (i.e., by replacing the 
delta function of Eq. (|70l) with its representation in the HO model space). These results 
are shown in Table [HI The model-space-restricted calculation is then similar, although 
not identical, to the standard NCSM calculation. In particular, as in the NCSM one 
loses the correct asymptotic behavior of the n-^^Be wave function. We observe that in 
the full NCSM/RGM calculation both relative kinetic and potential energies are smaller 
in absolute value. This is an effect of the re-scaling of the relative wave function in the 
internal region, when the Whittaker tail is recovered. The difference is significantly more 
substantial for the relative kinetic energy than for the potential energy. As a result one 
obtains a dramatic decrease of the energy of the 1/2+ state, which makes it bound 
and even leads to a g.s. parity inversion. This study shows that a proper treatment of 
the coupling to the n-^^Be continuum is essential in explaining the g.s. parity inversion. 
However, we cannot exclude that the NNN force plays a role in the inversion mechanism, 
until accurate calculations with both the NNN force and full treatment of the n-^°Be 
tail have been performed. 

8.4- Outlook: Ah initio no- core shell model with continuum 

It is possible and desirable to extend the binary-cluster [A— a, a) NCSM/RGM basis 
by the standard A-nucleon NCSM basis to unify the original ah initio NCSM and 
NCSM/RGM approaches. This will lead to a much faster convergence of the many- 
body calculations compared to the original approaches and, most importantly, to an 
optimal and balanced unified description of both bound and unbound states. 

In particular, we can generalize the expansion of the many-body wave function 
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given in Eq. (IM]) by explicitly including a set of A-nucleon NCSM eigenstates: 

I^J^T^ = ^cf ^lAAJ'^T) + drr'^^^A, IC^) , (73) 

where H^^^^\AXJ^T) = Ex\AXJ^T) with i/^^SM g^^^^ g_g_^ projected 
on the A^max^^ space or by the NCSM effective Hamiltonians ( |2^ or ( l29l) (with 
the -ffcM subtracted) also defined on the N^i,^hQ space. By projecting the many- 
body Schrodinger equation on the binary-cluster channel states ( l59l) and the NCSM 
eigenstates, we arrive at a system of coupled equations that can be schematically written 
as 

^NCSM 





Here, Ti. and A/" are the Hamiltonian and norm integration kernels defined in Eqs. ([H 
and ( l65i) . respectively. The q are the overlap functions introduced in Eq. ( l53i) and the 
f) are "vertex" functions defined by matrix elements {AXJ^T\HAu\^'l^.'^) with H the 
intrinsic Hamiltonian that can be expressed, e.g., as in Eq. (1661) . It is straightforward 
to implement this new approach that we name ab initio NCSM with the continuum 
(NCSMC). 



9. Conclusions 



The ab initio NCSM has evolved into a powerful many-body technique. In this review, 
we presented some of recent results obtained within this approach. We discussed, in 
particular, calculations with chiral EFT NN and NNN interactions for both s-shell and 
p-shell nuclei used, on the one hand, as a tool to determine the NNN interaction low- 
energy constants and, on the other hand, to predict properties of light nuclei. These 
calculations demonstrate the importance of the NNN interaction for nuclear structure. 
Recent advances in experimental techniques that allowed precise measurements of radii 
and moments of exotic isotopes motivated us to perform large basis NCSM calculations 
for He, Li and Be isotopes. These calculations were overviewed in this paper. We 
also described efforts to extend the NCSM calculations to larger model spaces and 
heavier nuclei by means of the importance-truncated calculations and by development 
of effective interactions for model spaces with a closed core. The most significant new 
developments, at least in our view, were discussed in the last part of this review. 
Extension of the NCSM to describe scattering and nuclear reactions via the RGM 
technique serves as a bridge to a development of a unified ab initio description of light 
nuclei with both bound and unbound states described simultaneously and treated on 
the same footing. Extensions of the NCSM/RGM formalism to include two-nucleon 
(deuteron), three- nucleon (triton and ^He) and four-nucleon (^He) projectiles are now 
under way. As a large HO basis expansion is needed in this formalism, not just for 
the convergence of the target and projectile eigenstates, but also for the convergence of 
the localized parts of the integration kernels, a combination of this approach with the 
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importance-truncated NCSM is key in a successful application of the NCSM/RGM and 
the NCSM with the continuum to heavier nuclei. 

There are other developments in the NCSM calculations that were not covered by 
this review that deserves attention. One of them is the recent analysis of NCSM wave 
functions by means of the representations of the symplectic Sp(3,R) group [1491 1150] 
and attempts to develop an NCSM code for calculations within SU(3) C Sp(3,R) 
symmetry-adapted basis. This approach aims to augment the model space by nuclear 
collective correlations, which are required for a microscopic description of monopole 
and quadrupole vibrational and rotational dynamics. This would allow, in particular, a 
reahstic description of giant quadrupole resonances. 

Further, the NCSM calculations were also used as a tool for development of the 
similarity renormalization group evolved NN interactions \118\ 1151] and the UCOM 
interactions [33j. Also, the NCSM method was recently adapted for systems of strongly 
interacting bosons in a trap |152j . 

Finally, there are remarkable advances in the development of the code MFD |81j . 
Calculations on tens of thousands of processors are now possible |153j and results 
in model spaces up to A^max = 10 for ^^C and ^^N with NN interactions are within 
reach |154] . Furthermore, calculations in the A^max = 8 model space with the NN and 
NNN interactions from Section [3] have already been performed for A = 7 and A = 8 
nuclei. 
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